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SUMMARY 


An engineering prediction method and associated computer 
code NOZVTX to predict nose vortex shedding from circular and 
noncircular bodies in supersonic flow at angles of attack 
and roll are presented. The body is represented by either 
a supersonic panel method for noncircular cross sections or 
line sources and doublets for circular cross sections, and 
the lee side vortex wake is modeled by discrete vortices in 
crossflow planes. The three-dimensional steady flow problem 
is reduced to a two-dimensional, unsteady, separated flow 
problem for solution. Comparison of measured and predicted 
surface pressure distributions, flow field surveys, and 
aerodynamic characteristics are presented for bodies with 
circular and noncircular cross sectional shapes. 

INTRODUCTION 


Current missile applications requiring high aerodynamic 
performance can involve noncircular body shapes in supersonic 
flow at high angles of attack and nonzero roll angles. The 
angle of attack range may be sufficiently high to cause 
formation of body separation vortices, and the body vortex 
shedding characteristics are directly influenced by the body 
cross-sectional shape and the flow conditions. It is desir- 
able to model the body vortex wake by means of a rational 
method capable of considering a variety of body shapes over a 
wide range of flow conditions. It is important that the 
separation vortex wake induced effects on the nonlinear 
aerodynamic characteristics of the body be handled properly. 

The phenomena of interest are the sheets of vorticity 
formed on the lee side of the body at high angles of attack. 
The vorticity is formed by boundary-layer fluid leaving the 



body surface from separation points on both sides of the 
body (Fig. 1) and rolling up into a symmetrical vortex pair. 

A method to predict these flow phenomena in the vicinity of 
circular and noncircular bodies in subsonic flow is described 
in References 1 and 2. The extension of the method to 
supersonic flow is described in Reference 3. 

The purpose of this report is to describe an engineering 
prediction method and associated computer code developed to 
calculate the nonlinear aerodynamic characteristics and flow 
fields of noncircular noses at high angles of attack in 
supersonic flow. The objectives of the method are to use a 
three-dimensional, attached flow, supersonic panel method, or 
a supersonic line source and doublet method to represent the 
body and a two-dimensional, incompressible, separated flow 
model to calculate the lee side vortex shedding from the body 
alone at angle of attack and angle of roll. The predicted 
pressure distribution on the body under the influence of the 
free stream and the separation vortex wake will be used to 
calculate the aerodynamic loads on the body. Conformal 
mapping techniques are used to transform noncircular cross 
sections to a circle for calculation purposes. 

The following sections of this report include a discus- 
sion of the approach to the problem and a description of the 
analysis and flow models required to carry out the calculation. 
The prediction method is evaluated through comparison of 
measured and predicted results on a variety of body shapes, 
including circular and elliptical cross sections. A user's 
manual for the computer code is also included. The manual 
consists of descriptions of input and output, and sample 
cases . 
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GENERAL APPROACH 


Bodies at high angles of attack exhibit distributed 
vorticity fields on their lee side due to boundary-layer fluid 
leaving the body surface at separation lines. One approach 
for modeling these distributed vorticity fields has involved 
the use of clouds of discrete potential vortices. Underlying 
the basic approach is the analogy between two-dimensional 
unsteady flow past a body and the steady three-dimensional 
flow past an inclined body. In fact, the three-dimensional 
steady flow problem is reduced to the two-dimensional, 
unsteady, separated flow problem for solution. Linear theory 
for the attached flow model and slender body theory to repre- 
sent the interactions of the vortices are combined to 
produce a nonlinear prediction method. The details of the 
application of this approach to the prediction of subsonic 
flow about circular and noncircular bodies are presented in 
References 1 and 2, and the extension to supersonic flow is 
described in Reference 3. Other investigators have also used 
this approach to successfully model the subsonic flow 
phenomena in the vicinity of circular cross section bodies 
(Refs. 4 and 5). The purpose of this report is to document 
the extension of the previous subsonic analysis of References 1 
and 2 to supersonic flow. 

The calculation procedure is carried out in the following 
manner. The body is represented by supersonic source panels 
(noncircular cross section) or supersonic line sources and 
doublets (circular cross section) , and the strength of the 
individual singularities is determined to satisfy a flow 
tangency condition on the body in a nonseparated uniform flow 
at ang^Le of attack and roll. Starting at a crossflow plane 
near the body nose, the pressure distribution on the body is 
computed using the full Bernoulli equation. The boundary 
layer in the crossflow plane is examined for separation 
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using modified versions of Stratford's separation criteria. 

The Stratford separation criteria are based on two-dimensional 
incompressible flow. At the predicted separation points, 
incompressible vortices with their strength determined by 
the vorticity transport in the boundary layer are shed into 
the flow field. The trajectories of these free vortices 
between this crossflow plane and the next plane downstream are 
calculated by integration of the equations of motion of each 
vortex, including the influence of the free stream, the body, 
and other vortices. The pressure and trajectory calculations 
are carried out by mapping the noncircular cross section 
shape to a circle using either analytical or numerical 
conformal transformations. The vortex-induced velocity con- 
tribution to the body tangency boundary condition includes 
image vortices in the circle plane. At the next downstream 
crossflow plane, new vortices are shed, adding to the vortex 
cloud representing the wake on the lee side of the body. 

This procedure is carried out in a stepwise fashion over the 
entire length of the body. The details of the individual 
methods combined into the prediction method are described 
in the following section. 


METHODS OF ANALYSIS 


The development of an engineering method to predict the 
pressure distributions on arbitrary missile bodies in super- 
sonic flow at high angles of incidence requires the joining 
of several individual prediction methods. In this section, 
the individual methods are described briefly, and the section 
concludes with a description of the complete calculation 
procedure . 
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Conformal Mapping 


The crossflow plane approach applied to arbitrary 
missile bodies results in a noncircular cross section shape 
in the presnece of a uniform crossflow velocity and free 
vortices in each plane normal to the body axis. The pro- 
cedure used to handle the noncircular shapes is to determine 
a conformal transformation for mapping every point on or 
outside the arbitrary body to a corresponding point on or 
outside a circular body. The two-dimensional potential flow 
solution around a circular shape in the presence of a 
uniform flow and external vortices is well known and has been 
documented numerous places in the literature (Refs. 4 and 6) . 
Thus, the procedure is to obtain the potential solution for 
the circular body and transform it to the noncircular body. 
Conformal transformations used are of two distinct types, 
analytical and numerical . 


Analytic transformation .- For very simple shapes like an 
ellipse [Fig. 2(a)], the transformation to the circle can be 
carried out analytically as described in Reference 7. For 
example , 


a 


v 



( 1 ) 


where 


a = y + iz 


( 2 ) 


in the real plane, and 


v = t + iX 


(3) 


in the circle plane. The derivative of the transformation 
is 


da 

dv 


1 



(4) 
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which is required for the velocity transformation discussed 
later. 

Numerical transformation .- For complex noncircular 
shapes, the transformation cannot be carried out analytically 
and a numerical transformation is required. The numerical 
transformation chosen is described in detail in Reference 8, 
and a brief summary of the conformal mapping procedure 
follows . 

The sequence of events in the numerical mapping is shown 
in Figure 2 (b) . The arbitrary cross section shape of the 
body in the a-plane is required to have a vertical plane of 
symmetry. The transformation of interest will map the 
region on and outside the body in the a-plane to the region 
on and outside a circle in the v-plane. The first step is 
a rotation to the C -plane so that the cross section is 
symmetric about the real axis. Thus, 

o ~ y + iz (5) 


and 


C = io 

= K + in (6) 

A mapping that transforms the outside of the body in the 
C-plane to the outside of the unit circle is 





m 


TT<v 

r=l 



dv 


(7) 


where y r are the exterior angles of the m segments of the 
body cross section. In Equation (7) , tt denotes a product 
series. For a closed body 
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( 8 ) 


m 

I Y r = 2 tt 
r=l r 

The final transformation has the form 


a 


-l 



+ 


L 


I 

1=1 



(9) 


where the A^ coefficients are obtained through an iterative 
scheme described in Reference 8 and r Q is the radius of the 
equivalent circle in the v-plane. The derivative of the 
transformation is 


da 

dv 



L 


I 

1=1 


r ^ +1 lA 
r o 

Z+I 

v 


( 10 ) 


which is required for the velocity calculations described in 
a later section. 

The above numerical mapping procedure has been applied 
to a wide range of general cross section shapes with good 
success . 


Body Model 

The prediction method described in this report is 
applicable to both circular and noncircular bodies. Separate 
body models are developed for circular and noncircular shapes 
because arbitrary noncircular shapes require panels and 
numerical conformal mappings and circular bodies can be 
represented by more computationally efficient means. Both 
flow models are described in this section. 

Circular bodies .- A potential flow model of an axisym- 
metric body in supersonic flow is described in References 9 
and 10. Such a body is represented by a distribution of line 
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sources and sinks and line doublets on the body axis to 
account for volume and angle of incidence effects, respec- 
tively. The strengths of the line singularities are 
determined from the flow tangency condition applied at points 
on the body surface. A downstream marching procedure from 
the nose tip to the body base is used for this calculation. 

The axis singularity distribution for bodies of revolu- 
tion provides an accurate and efficient means to calculate 
the induced velocity field associated with the body at 
angle of incidence. The calculation requires much less 
computation time than a panel method, and the resulting 
velocity distributions are quite smooth. One restriction on 
the body shape is the requirement that the nose contour lie 
inside the Mach cone. This constitutes the major limitation 
to the method . 

Noncircular bodies .- The configurations considered in 
this study are bodies alone, without fins or control surfaces. 
The body geometry is described by the coordinates of the 
cross-sectional shape at a number of axial stations, and 
linear interpolation between specific points is used. The 
body is represented by a three-dimensional lattice of super- 
sonic source panels distributed on the surface of the body. 

The panels account simultaneously for volume and angle of 
incidence effects. The attached flow solution for a panel 
inclined to the flow direction is based on supersonic linear 
theory as described in Reference 11 and the method used is a 
modified version of the program code presented in 
Reference 12. The pitch and yaw angles are related to the 
included angle of attack and angle of roll by the expressions 

sina = sina cos<b {11) 

c 

sing = sina c sin<|) (12) 
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The flow tangency condition is applied at the control 
point located at the panel centroid of each body panel . The 
result is a set of linear simultaneous equations in terms of 
the unknown panel source strengths. Solution by an iterative 
approach produces the singularity strength on each panel and 
thus provides a means to calculate the velocity field induced 
by the body. The nature of the solution is such that the 
velocity on the body surface is correct only at the control 
points; therefore, linear interpolation is used to determine 
body surface velocity components for use in pressure calcula- 
tions at points lying between control points. The body panel 
induced velocities are calculated one time for the entire 
body to provide a table of velocity components for interpola- 
tion purposes. Using a slender-body theory approximation, 
these velocity components are combined with other velocity 
components for the pressure calculation described in a later 
section . 

There is a basic limitation in the use of the above 
panel method. The angle of inclination between a body source 
panel and the body centerline must be less than the semi- 
apex angle of the Mach cone associated with the free-stream 
Mach number. Thus, there is a limit to the body nose con- 
figurations which can be modeled by the body source panel 
method . 


Vortex Shedding Model 

The body vortex shedding model described in this section 
is nearly the same as the subsonic model presented in 
Reference 2; therefore, the following description will be 
brief, and only the differences between the two models will 
be discussed in detail. 
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Equations of motion .- The equations of motion of a shed 
nose vortex in the presence of other free vortices in the 
vicinity of a body in a uniform stream follow. In the circle 
(v) plane, the position of a vortex, r , is 



t + iA,„ 
m m 


and the image of F m is located at 



In the real plane, the position of the vortex r is 


r = y + iz 
J m m 


(13) 


(14) 


(15) 


The complex potential in the real plane is 


W (a) = $ + iY 


(16) 


and the corresponding velocity at is 


dw a) , 

m d . . , dv 

m m da dv m da 


(17) 


o=o 

i 

v=v 


m 

) 

m 


The complex potential of T is not included in Equation (17) 
to avoid the singularity at that point. The derivative of 
the transformation is obtained from Equation (4) or 
Equation (10) . 

The total velocity at r in the crossflow plane is 
written as 


m 


- iw, 


m 


V 


G + G a + G + G + G_ + G 
a 3 n m T r 


(18) 


10 



where each term in Equation (18) represents a specific 
velocity component in the a-plane. The first term represents 
the uniform flow due to angle of attack. 


G 


a 




r ^ 

2 “ 


i sina 

1 + 

r 

o 

V 


dv 

da 



l m J 




a=o 


m 


(19) 


The second term represents the uniform flow due to angle of 
yaw . 


G e 


-sin$ 



r ^ 

r 

o 

2" 

dv 

1 - 

V 

m 


da 


k, j 

. 



( 20 ) 


The next term represents the influence of all vortices and 
their images, with the exception of r . 


N r 

_ _ . y n 

n - 1 2iTr V 

n=l o 00 


( V r o ) " ( VV 


(v m /r o>- (v n /r o> 


dv 

da 


n^m 


( 21 ) 


a=a 


m 


The next term is due to the image of f m . 


r m r 

_ ^ « in 

G m = 1 27rr V 


O oo 


L (v m /r o ) - (r o /v m ) J 


dv 

da 


( 22 ) 


o=o. 


m 


The fifth term in Equation (18) represents the potential 

of r in the a-plane and is written as 
m 


g t = - 1 


m 

2ttV 


'll d 

dvl 

L 2 J dv 

|d°J 


0=0 


m 


(23) 
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The last term in Equation (18) represents the velocity com- 
ponents induced by the portion of the body panel singularities 
representing the body volume effects. 


G 


r 



V 

00 


(24) 


Since the body panel singularities are three dimensional, 
they contribute an induced axial velocity, u . 

The differential equations of motion for are 



da 

v - iw 


m 

m m 


dx V ( 

cosa + u 
» c r 

where 

a 

= y - iz 


m 

J m m 

Therefore, the two 

equations 

which must be integrated 

the body length to 

determine 

the trajectory of are 


dy 
J m 

V 

m 


dx V 

< 

cosa + u 
» c r 

and 

dz 

w 


m 

m 


(25) 


(26) 


dx V cosa + u 
» c r 


(27) 


(28) 


There are a pair of equations like (27) and (28) for each 
vortex in the field. As new vortices are shed, the total 
number of equations to solve increases by two for each added 
vortex. These differential equations are solved numerically 
using a method which automatically adusts the step size to 
provide the specified accuracy. 


Surface pressure distribution .- The surface pressure dis- 
tribution on the body is required to calculate the forces on 
the body and the separation points. The surface pressure coef- 
ficient is determined from the Bernoulli equation in the form 
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Y 

Y-l 


(29) 


where 



2 




+ 


Y-l 

2 


(C, 


and 




U 


V 


2cosa 


V. 


dx 


(30) 


(31) 


In Equation (31) , U is the total velocity (including V^) at 
a point on the body. It consists of many of the components 
described in the previous section. The body induced velocity 
components from Equation (24) include the total effect of the 
panel solution at angle of attack and roll. As a consequence 
of using the full panel solution in this component, the two- 
dimensional doublet part of Equations (19) and (20) are not 
required. The vortex wake contribution is included through 
Equations (21) and (22) . 

The last term in Equation (31) represents the axial 
velocity missing from the two-dimensional singularities 
in the flow model. In this case, the shed vortices are the 
only singularities contributing to this term. The complex 
potential of a number of vortices, r , and their images in 
the circle plane is 


W ( v ) = <f> + i ¥ 


N 


= I 

n+1 



r 2i 



r 


-i £n(v - v ) + i £n 

v - -2- 


V 



n J 



(32) 
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The "unsteady" term in Equation (31) , evaluated on the 
body surface, is 


= Real sm. 

dx dx 


r=r 


(33) 


Equation (33) becomes 


d 4> 
dx 


N 

I 

n=l 


r 

n 

2tt 


** 

J 

dX 

(x - x ) , n - 

n dx 

(X - 

V 

dx 

n 

dx 

\ 

(x - x ) 2 + 
n' 

(X - 

V 

2 


, 2 2 N 

ir — x r ) 
n no' 


2Ax 


dx dA 0 dX dr 

+ 2XX -=~2 - r 2 — ~ - 2r X -=-£ 
n dx n dx o dx o n dx 


2 2 2 2 2 2 
(xr - x r ) + ( Ar - X r^) 

n no n no 


+ 


2 2 
< Xr n' Vo> 


dX dx 

2x A + 2xx -=-2- 

n dx n dx 


- dx dr 

2 n ~ o 

- — “ 2r x -j— i 

o dx o n dx J 


2 2 2 2 2 2 
(xr - x_r ) + (Ar - A r ) 

n no n no 


(34) 


where 


2 2,2 

r = x + A 
n n n 


35) 


Separated wake .- The separated wake on the lee side of 
the body is made up of a large number of discrete vortices, 
each vortex originating from one of two possible locations 
at each time step in the calculation. The major portion of 
the lee side vortex wake has its origin at the primary 
separation points on each side of the body. The remainder 
of the wake originates from the secondary separation points 
located in the reverse flow region on the lee side of the 
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body. Both of these points are illustrated in the sketch of 
a typical crossflow plane of an elliptic cross section body 
shown in Figure 3. The mechanics of the calculation of the 
individual vortices follows. 

As described in References 1 and 2, the pressure dis- 
tribution for the primary flow in the crossflow plane is 
referenced to the conditions at the minimum pressure point, 
and a virtual origin for the beginning of the boundary layer 
is assumed. The adverse pressure distribution downstream of 
the minimum pressure point is considered with either 
Stratford’s laminar (Ref. 13) or turbulent (Ref. 14) separa- 
tion criterion to determine whether or not separation has 
occurred. These criteria, based on two-dimensional, 
incompressible, flat plate data, are adjusted for three- 
dimensional crossflow effects in Reference 2. The laminar 
separation criterion states that the laminar boundary layer 
separates when 

- 0.087 sina c (36) 


, 1/2 


_E 


In a turbulent boundary layer, separation occurs when 

, 1/2 


dC 




-6 - 0 * 1 

(Re^ x 10 ) - 0.35 sina ( 


(37) 


If the criteria indicate a separation point, the 
vorticity flux across the boundary layer at separation is 
shed into a single point vortex whose strength is 




Ax 

cosa 

c 


(38) 


assuming no slip at the wall. 
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The initial position of the shed vortex is determined 
such that the surface velocity at the separation point is 
exactly canceled by the shed vortex and its image. When this 
criterion results in a vortex initial position that is too 
near to the body surface, certain numerical problems can 
cause difficulty in calculating the trajectory of this 
vortex. If the intial position of the vortex is nearer than 
five percent of the body radius from the surface, the vortex 
is placed five percent of the equivalent body radius, from the 
surface. This position represents the approximate thickness 
of the boundary layer. 

The calculation of secondary separation is carried out 
in the same manner as described above for primary separation. 
It is necessary that a reverse flow region exist on the lee 
side of the body and that a second minimum pressure point be 
found in this region. Surface flow visualization of 
secondary separation regions on bodies at high angles of 
attack in subsonic flow are shown in Reference 15. For 
purposes of this analysis, the reverse flow is assumed to be 
laminar from the lee side stagnation point to the secondary 
separation point, and Stratford’s laminar criterion is used 
to locate the secondary separation point. Laminar separation 
in the reverse flow region is expected because of the low 
velocities on this portion of the body. The vortex released 
into the flow at the secondary separation point has the 
opposite sign of the primary vortex and is generally much 
weaker in strength. 

Forces and moments .- The forces and moments on the body 
are computed by integration of the pressure distribution 
around the body. At a specified station on the body, the 
normal-force coefficient on a Ax length of the body is 
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2tt 

r 


C^r' cosB' d3 


(39) 


c 

n 


'An' 
q D 

■*00 


1 

D 


where r 1 is the distance from the axis of the body to the 
body surface and 3 1 is the local slope of the body in the 
crossflow plane. This is illustrated in the sketch in 
Figure 4. For circular bodies, r f = and 3 1 =3. The 
total normal force coefficient on the body is 


C 


N 




c n dx 


(40) 


and the pitching-moment coefficient is 


M 


C = 


m q S 




ref 


L 

'x - x' 
m 

J n 

o 

~Z 7“ 
ref J 


dx 


(41) 


The center of pressure of the normal force is 


x 


cp 


'ref 


n m 

I 


ref 


m 


'N 


(42) 


Similarly, the side-force coefficient on a Ax length 
of the body is 


c = 


As 
( Ax 


2tt 


q D 

^■oo 


1 

" " D J 


C r 1 sin3 1 dB 
P 


The total side-force coefficient on the body is 


c y q„s s j 


c dx 

y 


(43) 


(44) 
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and the yawing -moment coefficient is 


M 


r — r 

^n cTsZ 7 

oo ref 


±j 

•-u 


X_ - xl 
m 


'ref 


dx 


The center of pressure of the side force is 


(45) 


x 

I 


cp C 

3: = _n + 

ref y 


x 


m 


ref 


(46) 


A noncircular body at an arbitrary roll angle may 
experience a rolling moment caused by the nonsymmetry of the 
loading around the body. The total rolling moment on the 
body is calculated by summing the moments of the individual 
components of normal force and side force around each cross 
section and integrating over the body length. Using 
Equations (39) and (43) , the rolling moment coefficient is 
calculated as 


M 


C l q SI 


x 


ref 


2t t 
r 


SI 


ref 


(yC^r f cos$ 1 


zC r' sin3')d$|dx 
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Vortex smoothing .- A problem common to discrete vortex 
models is the singularity in the vortex induced velocity 
distribution at the loction of a point vortex and the very 
large velocities induced near the vortex. This can lead 
to errors in the surface pressure distribution when there 
are vortices near the body surface and to errors in the 
vortex tracking when vortices pass in close proximity. 
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The vortex diffusion model described in Reference 2 has the 
effect of removing the singularity at the vortex, but very 
large induced velocities are still in evidence near the 
vortex, thus the noted roughness in the surface pressure 
distribution continues to be a problem. Several solutions to 
this problem have been investigated, and a discussion of 
these vortex flow models follows. 

In Reference 16, discrete vortices are replaced by a 
distribution of vortices on straight line segments connect- 
ing the discrete vortex positions. In this way, the 
algebraic singularity is reduced to a logarithmic singularity. 
The strength of the uniform distribution of vorticity on each 
line segment is determined from the strengths of the discrete 
vortices at the ends of the segment. Defining as the 
complex velocity induced at v m by the line segment, analogous 
to Equation (21) , the velocity becomes 


’d. 

3 


ir n' 

V n+1 

- V 

n 

log 

r \ 

- V 

n+1 m 

27TV 

oo 

1 V 

- v 1 

v - V 

n m 


1 n+1 

n 1 

> 


v ./ 


(48) 


where v m is not on the line segment. The strength of the 
vortex segment is 
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except when T n or represents the last vortex in the 

chain. When this occurs, the remainder of the strength of 
the last vortex is included in the last line segment to 
avoid losing a portion of the shed vorticity. 

This smoothing procedure was investigated for a number 
of calculations, and it accomplished the smoothing of the 
induced velocity field to a great extent. The surface 
pressure distributions had some residual roughness caused 
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by other effects, but the pressure spikes caused by locally 
high vortex-induced velocities were eliminated. The major 
disadvantage of the above smoothing technique is cost. 

The nature of the total calculation procedure is such that 
many velocity calculations are required; therefore, the 
additional effort to set up the smoothing procedure and 
calculate velocities at numerous body points causes a 
dramatic increase in computation time. For this reason, 
alternate smoothing techniques are considered. 

The smoothing technique incorporated into the present 
version of NOZVTX is one in which a single point vortex 
is represented by a simple distribution of several point 
vortices. The smoothing calculation is carried out only 
for vortices which lie within 10% of the local body radius 
of the surface of the body, and it is used only for calcu- 
lating pressures on the body surface. The vorticity is 
distributed according to the procedure shown in Figures 5 (a) 
and (b) . In Figure 5 (a) , vortex 2 is the only vortex which 
lies within 10% of the local body radius of the surface. 

The "distributed vorticity" model consists of nine smaller 
vortices; one located at the actual position of vortex 2, 
four located along the line joining vortices 1 and 2, and 
four along the line joining vortices 2 and 3. Each vortex 
has a strength equal to one-ninth of the actual vortex 
strength. Figure 5(b) shows distributed vorticity models 
for vortices which are end points of the vortex cloud; 
that is, for the first vortex and the most recently shed 
vortex. In the case of the first shed vortex, r^, the 
distributed vorticity model consists of five point vortices, 
each vortex having a strength equal to one-fifth that of the 
actual vortex. The model for the most recent shed vortex, 
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r 3 , is shown in Figure 5(b). The vorticity is distributed 
along two separate lines. One line is defined by the 
actual vortex position and its shed position at the 
previous axial station, indicated by 3* in Figure 5(b) . 

The second line is defined as before, the end points being 
the locations of T 2 and r^. 

The smoothing technique described above was investigated 
for a number of geometry and flow condition combinations, and 
it was found to smooth pressure spikes caused by locally 
high vortex-induced velocities or values of d^/dx. As 
expected, use of the smoothing technique increases execution 
time, but the actual execution time is determined by the 
number of vortices which lie close to the body surface. 

This technique redistributes the vortices and tends to spread 
their effect over a larger area, but the model is still 
made up of discrete vortices, each of which has a singularity 
at its origin. A viscous core model for the individual 
vortices is required to remove the singularity in the induced 
velocity calculation. This model is described in the 
following section. 

Vortex core .- The diffusion core model (Ref. 2) for the 
point-vortex induced velocities removes the singularity at 
the vortex origin and effectively reduces the velocities near 
the vortex. The tangential velocity induced by a single 
point vortex is written as 
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where r is the distance from the vortex to the field point 
and x is a measure of the age of the vortex. The induced 
velocity from this core model is illustrated in the 
following sketch. 



Vortex Induced Velocity 


The vortex core model represented by Equation (50) 
has received considerable attention in the context of body 
vortex induced effects, and it has a number of shortcomings. 
Since the exponential term is a function of r, the flow 
medium (v) , and the age of the vortex (x) , the core is 
constantly changing size as the vortex moves through the 
field. Under certain conditions, the core radius, denoted 
as r c in the sketch, can become very small and the induced 
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velocity becomes unrealistically large. In an attempt to 
further modify the core model to keep the induced velocities 
within physically realistic limits, the following modifica- 
tion was made. The location of the maximum induced velocity, 
r c in the sketch, is fixed at a specific radius to be 
selected by the user. Given a core radius, the vortex 
induced velocity is 




f 2 


1 - exp 

-1.256 



l cj 



(51) 


The core model described by Equation (51) is included 
in the version of the code described in this report, and in 
a later section, guidelines for the selection of an appro- 
priate core radius are presented. Results indicate that this 
simple core model provides adaquate smoothing for the 
necessary velocity calculations. 

Other investigators of discrete vortex models must use 
a core model of some type. Although there are many different 
core models, all serve the same purpose as those described 
above, and nearly all are directed at eliminating the 
singularity and reducing the maximum induced velocity. The 
code described in a later section is easily modified to 
incorporate another core model if the user desires to do so. 


PROGRAM 

General Description 

The computer code, NOZVTX, presented in this section 
has application as an engineering prediction method for 
circular and noncircular bodies at high angles of attack 
in supersonic flow. The code is in effect a combination of 
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two separate codes. The first is the subsonic body vortex 
shedding method described in References 1 and 2, and the 
second is the supersonic body model for bodies alone described 
in References 10, 11, and 12. The combined code provides 
a means to predict the nonlinear aerodynamic characteristics 
and flow fields for bodies with circular and noncircular 
cross sections. The effect of fins is not included. 

The code consists of a main program, NOZVTX, and 50 
subroutines. The overall flow map of the program is shown 
in Figure 6 where the general relationship between the sub- 
routines and external references can be seen. An alpha- 
betical listing of the subroutines is shown in Table I . 

It was an objective of this effort to keep the two codes as 
distinct as possible to maintain their individual identity 
so that modifications to either code could be carried out 
without major changes in the other code. Note that the 
main subroutine of the body vortex shedding calculation, 

ELNOSE , includes the body panel code via the calls to 
subroutines GEOM , VELCMP , and SOLVE. The line singularity 
model for axisymmetric bodies is included through a call to 
subroutine BDYGEN. The communication of information between 
segments of the program is handled by name common blocks. 

A cross reference table for the calling relationship 
between the program subroutines and the external references 
is shown in Figure 7. A similar table for the common blocks 
is shown in Figure 8. 

The program is written in standard FORTRAN language. 
Operation requires four (4) auxiliary files in addition 
to the usual input and output units. Core storage for 
execution is approximately 230,000 octal on a CDC CYBER 760 
computer. Execution times on the same computer can vary 
from 10 to 400 or more seconds. Execution time depends on 
many factors such as the number of panels required for the 
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body, the integration interval (DX) , the number of shed 
vortices, the use of symmetry, and even the flow conditions. 

Limitations 

Program NOZVTX is applicable to bodies alone at high 
angles of attack in supersonic flow, with some limitations. 
Certain of the limitations to follow are basically rules-of- 
thumb for the use of the program that have evolved through 
experience with the program and verification of the program 
by comparison with data. The range of applicability of the 
program may be extended by additional verification with 
appropriate data. 

There is a Mach number /geometry limitation in the 
body model that requires the body surface to lie behind the 
Mach cone. This has greatest effect at high Mach numbers 
or on blunt noses at low Mach numbers, and a message is 
printed if the geometry condition is violated. It has been 
the experience of the authors that, when this condition 
arises, the nose can generally be made slightly more pointed 
so that it lies within the Mach cone. Often this modifica- 
tion is required only on the first few percent of body 
length, and the effect on the total nose forces and moments 
and shed vortex field is minimal. 

An associated Mach number limitation involves the 
quality of the predicted surface pressure distribution at 
high Mach numbers. As will be shown in comparisons of experi- 
ment and theory, the predicted pressure coefficients on the 
windward side are much lower than those measured at body 
points which are near to the shock wave. This is a common 
problem for linear body models like the axis singularities 
and the panels considered in NOZVTX, and it is possible that 
the windward pressure can be corrected by use of nonlinear 
pressure prediction techniques. 
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Another less rigid limitation involves the Mach number/ 
angle of attack range. It is suggested that the crossflow 
Mach number (M sina ) not exceed unity. A stricter limita- 
tion is that the crossflow Mach number should not exceed the 
critical value; that is, the local Mach number in the cross- 
flow plane should not exceed unity at any point in the plane. 
The root of this limit is the vortex shedding calculation 
which uses an incompressible vortex model and an incompress- 
ible separation criterion. Also, when the crossflow Mach 
number exceeds the critical value, a complex shock near the 
shoulder may contribute to separation. It is important that 
the user be aware of these limits, particularly when the 
program is being used in that range. However, the author has 
used the program for flow conditions in which the crossflow 
Mach number slightly exceeded unity, and the results were 
quite resonable. Some of these results are presented in a 
later section. 


Error Messages and STOPS 

The program code has numerous internal error messages. 
These are generally explicit and they will be described 
below only if an explanation will provide some insight into 
the problem and its correction. There are several execution 
terminations at numbered STOP locations within the program. 
These are described in the table below along with some 
suggestions to correct the problem. 

The most common message is "CONVERGENCE NOT OBTAINED IN 
INTEGRATION." This message, along with the axial station, 

XI, and the integration interval, H, is printed when the 
trajectory calculation does not converge within the specified 
tolerance, E5. One message will terminate execution at 
STOP 40 described in the table to follow. 
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The program contains several error messages pertaining 
to the number of panels used to describe the body. The 
message "ERROR - NUMBER OF ROWS OF PANELS IN BODY SECTION 1 
EXCEEDS 30" is self explanatory and results in STOP 20. A 
similar message "ERROR - NUMBER OF ROWS OF SINGULARITY 
PANELS ON BODY EXCEEDS 30" results in STOP 50. The message 
"ERROR - NUMBER OF BODY PANELS EXCEEDS 600" results in 
STOP 60. The message "ERROR - BODY HAS MORE THAT 20 
CIRCUMFERENTIAL PANELS" refers to the number of panels 
describing the body cross section circumference. This error 
results in STOP 200. All of the errors associated with the 
body panels are caused by incorrect values being input in 
items 32 and 41 of the input list. 

If a pivot element of the matrix of panel influence 
coefficients is zero, the message "ERROR - THE MATRIX IS 
SINGULAR" is printed and the program terminates execution 
at STOP 300. This is a very uncommon error and is not seen 
for normal body panel arrangements. 

As mentioned in a previous section, the slope of any 
panel on a noncircular body must be less than the Mach 
angle. If a panel has a greater slope, the message 
"ERROR - BODY PANEL SLOPE EXCEEDS MACH ANGLE" is printed 
and execution terminates at STOP 210. A modification of 
the body geometry to decrease the panel slopes is required 
to correct this situation. 

The surface of an axisymmetric body must lie inside 
the Mach cone from the nose tip. If certain body points lie 
outside the Mach cone, subroutine BDYGEN will make an effort 
to correct the definition of the body and revise the radius 
distribution. If the Mach number is too high to correct 
the body within its specified length, the message "BODY 
RADIUS TOO LARGE IN RELATION TO BODY LENGTH" is printed and 
execution terminates at a STOP. If, in the process of 
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correcting the body points outside the Mach cone, the 
subroutine cannot locate an appropriate point, the message 
"MACH CONE - BODY MERIDIAN INTERSECTION NOT FOUND AFTER 
100 TRIALS" is printed and execution terminates at a STOP. 
These difficulties are corrected by modifying the body 
shape or lowering the Mach number of the run. 

The numerical conformal mapping procedure can generate 
several error messages from subroutine CONFOR. The message 
" * * *ERROR IN SUM OF EXTERIOR ANGLES***" is printed if the 
routine has difficulty calculating appropriate body angles 
within a specified tolerance. Another message " ***COEFFICIENT 

AN (1 ) = GREATER THAN +- .001 AFTER ITERATIONS OUT 

OF " is printed on rare occasions when the mapping procedure 

is having difficulty converging on a solution. Both of the 
above messages are informational, and they do not result in 
a program STOP. However, the user is advised to check the 
transformed shape carefully with the actual shape before 
using the results. If there is a difficulty in the numerical 
mapping calculation, refer to the Input Preparation section 
for assistance. 


Numbered STOP Locations 


Termination 

Message 


Description 


STOP Normal program stop or the axisymmetric body 

lies inside the Mach cone from the nose tip. 

STOP 20 More than 30 rows of panels on body. Revise 

input values of NFORX and KFORX in input 
items 32 and 41, respectively. 

STOP 40 Vortex trajectory integration has failed to 

converge at an axial station. This problem is 
usually caused by (1) two or more vortices in 
close proximity such that they are rotating 
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about each other so rapidly that the integration 
subroutine cannot converge on a solution, or 
(2) a vortex too near the body surface. These 
situations are often dependent on the body shape 
and the specific flow conditions, but the follow- 
ing steps should be tried to resolve the 
problem. Increase the integration tolerance E5 , 
or increase the size of the combination param- 
eter RGAM so that it is greater than the distance 
between the troublesome vortices. It often 
helps to increase the integration interval DX 
and restart the calculation from a previous 
station . 

STOP 41 The vortex trajectory calculation has resulted 

in a vortex inside the body. This is a rare, 
but serious, problem; and it signifies a major 
problem with the solution. It usually means 
that the body paneling scheme is in error. If 
no obvious problem exists, change the paneling 
layout to increase the number of circumferential 
and axial panels. 

STOP 44 Total number of shed vortices exceed 200. 

Combine some vortices and restart the calcula- 
tion and run to completion. An alternate 
solution is to increase the integration interval 
DX so that fewer vortices are shed. 

STOP 45 The number of positive or negative primary 

separation vortices exceeds 70. Use the same 
remedy as for STOP 44. 

STOP 46 The number of positive or negative secondary 

or reverse flow separation vortices exceeds 30. 
Combine the secondary vortices to reduce their 
total number and restart the calculation; 
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STOP 50 
STOP 60 
STOP 200 

STOP 202 

STOP 210 
STOP 300 


increase DX to reduce numbers of vortices shed, 
or set NSEPR=0 and eliminate the secondary 
vortices from the calculation. 

More than 30 rows of panels on body. See 
STOP 20. 

More than 600 panels on body. Revise requested 
number of panels in input items 32 and 41. 

More than 20 panels on the body circumference. 
Revise values of NRADX and KRADX in input items 
32 and 41, respectively. 

There is an error in the vortex tables in sub- 
routine VTABLE . Check input values of NBLSEP , 
NSEPR, NVP , NVM, NVA , NVR, and RGAM in items 
1 and 6 . 

A body panel has a slope greater than the Mach 
angle . 

There is a singularity in the panel influence 
coefficient matrix. Check panel layout geometry 
for errors. 


An additional comment is required here about a situation 
which can occur which results in a fatal execution error and 
termination at a TIME LIMIT with no other message. This is 
not a usual occurrence, but it has been observed in certain 
special cases. The typical case is near the end of a very 
long body from which a large number of vortices have been 
shed. The difficulty is that two differential equations of 
motion must be integrated for each vortex and the derivatives 
involve all vortices in the field. Solution by a Kutta- 
Merson scheme can involve many derivative calculations . The 
minimum time required for each axial step in the calculation 
is proportional to the square of the total number of 
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vortices. As the number of vortices increases, convergence 
becomes more difficult and the time for each step can 
increase dramatically. If this situation should occur, the 
best approach is to combine adjacent vortices and restart 
the calculation with a larger convergence criterion E5. It 
is sometimes useful to increase the axial interval DX. 

Input Description 

The purpose of this section is to describe the input to 
program NOZVTX. The input formats are shown in Figure 9, and 
the definitions of the individual variables are provided in 
the following paragraphs . The input to the body panel and 
axis singularity portion of the program is maintained 
nearly identical to the input requirements for program WDYBDY 
of Reference 10. The purpose for this was to minimize 
changes in WDYBDY; however, it leads to some redundancy in the 
preparation. The small amount of redundancy is acceptable 
in light of the simplification in the use of the current 
program. If there is a question about the input in the panel 
portion of the program, the user is referred to Appendix K 
of Reference 10 for additional information. 

The remainder of this section includes a description of 
each of the input variables and indices required for the use 
of program NOZVTX. The order follows that shown in Figure 9 
where the input format for each item and the location of 
each variable on each card is presented. Data input in 
I-format are right justified in the fields, data input in 
F-format may be placed anywhere in the field and must include 
a decimal point, and data input in E-format are right justi- 
fied in the fields and must include the decimal point and 
the complete exponent designation. 
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Note that all length parameters in the input list are 
dimensional variables? therefore, special care must be taken 
that all lengths and areas are input in a consistent set of 
units . 

Item 1 is a single card containing sixteen integers, 
each right justified in a five-column field. The variables 
are defined as follows. 


Item 

Variable 


Description 

1 

NCIR 

Cross section 
=0, circular 
=1, elliptic 
=2, arbitrary 

shape index . 
cross section 
cross section 
cross section 


NCF 


ISYM 


NBLSEP 


NSEPR 


Numerical transformation index. 

=0, mapping coefficients are calculated 
based on body shape input 

=1 , mapping coefficients are input by 
user in items 18 through 22 

Symmetry index. 

=0, right-left flow symmetry 

=1, no symmetry (required if <J> ^ 0°) 

Body vortex separation index. 

=0, no separation (required if a c - 0° 
in item 5) 

=1, laminar separation 

- 2 , turbulent separation 

Reverse flow separation index. 

=0, no separation 

=1, laminar separation in reverse flow 
region 
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NS MOTH 


Vortex induced velocity smoothing index. 

=0 , no smoothing 

=1, smoothing of vortex-induced velocities 
used in pressure calculation 
=2, smoothing of vortex-induced velocities 
used in pressure calculation and 
smoothing of 3<j)/9x calculation 

NDFUS Vortex core model index. 

=0, potential vortex 

=1, diffusion core model; preferred (see 
RCORE in item 5) 

NDPHI Unsteady pressure term index. 

= 0, omit 34>/9x from Cp calculation 
=1, include 3 4>/3x term; preferred 

INP Nose force index. 

=0, slender body theory force on portion 

of nose ahead of starting point (XI) 
=1, zero force on nose ahead of XI; 

preferred value for restart option 

NXFV Number of x-stations at which flow field 

is calculated or special output generated. 
See item 23. 

0 £ NXFV <_ 8 

NFV Number of field points for flow field 

calculation. See item 24. 

0 < NFV <_ 200 

NVP Number of +T vortices on +y side of body 

to be input for restart calculation. See 
item 26. 

0 < NVP <70 
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NVR 


Number of -r reverse flow vortices on +y 
side of body to be input for restart 
calculation. See item 27. 

0 <_ NVR <_ 30 

NVM Number of -T vortices on -y side of body to 

be input for restart calculation. See 
item 28. NVM = 0 if ISYM = 0. 

0 £ NVM £ 70 

NVA Number of +T reverse flow vortices on -y 

side of body to be input for restart 
calculation. See item 29. NVA = 0 if 
ISYM = 0. 

0 <_ NVA £ 30 

NASYM Asymmetric vortex shedding index. See 

item 7 . 

=0, no forced asymmetry; preferred 
=1, forced asymmetry 

Item 2 is a single card defining seven integer output 
option indices. Each index is right justified in a five 
column field. The indices are defined as follows. 

Item Variable Description 

2 NHEAD Number of title cards in item 3 . 

NHEAD > 1 

NPRNTP Pressure distribution print index. 

= 0, no pressure output except at special 
x-stations specified by item 23 
=1, pressure distribution output at each 
x-station 
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NPRNTS 


NPRNTV 


NPLOTV 


NPLOTA 


NPRTVL 


Vortex separation print index. 

=0, no output 

=1, output at each x-station; preferred 

= 2, detailed separation calculation 

output; for debugging purposes only 

Vortex cloud summary output index. 

=0, no vortex cloud output 

=1 , vortex cloud output; preferred 

Vortex cloud printer-plot option. 

=0, no plot 

=1 , plot full cross section on a constant 
scale 

- 2 , plot upper half cross section on a 
constant scale 

=3, plot full cross section on a variable 
scale; preferred 

Plot frequency index. 

=0, no plots 

=1, plot vortex cloud at x-stations 
specified by item 23 

= 2 , plot vortex cloud at each x-station 

Velocity calculation auxiliary output for 

debugging purposes only (Sub. VLOCTY) . 

=0, no output; preferred 

=1, print velocity components at field 
points; see items 23 and 24 

— 2 , print velocity components at body 
control points during pressure 
calculations. This option can produce 
massive quantities of output (not 
recommended) . 
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Item 3 is a series of cards containing hollerith informa 
tion indentifying the run. The cards are reproduced in the 
output just as they are input. 


Item 

Variable 

Description 

3 

TITLE 

NHEAD cards of identification. Information 
to be printed at top of output sheets. 

Item 4 is a single card containing reference information 
used in forming aerodynamic coefficients. The variables are 
defined as follows. 

Item 

Variable 

Description 

4 

REFS 

Reference area. REFS > 0 


REFL 

Reference length. REFL > 0 


XM 

Moment center. 


SL 

Body length (L) . 


SD 

Body maximum diameter, or equivalent non- 
circular body diameter (D) . Used in 
definition of RE. 

Item 5 is a single card containing the flow condition 
parameters. Each variable is input in a ten column field. 
The flow parameters are: 

Item 

Variable 

Description 

5 

ALP H AC 

Angle of incidence, degrees. 

0° < a < 90° 

— c 

If a =0°, NBLSEP = 0 in item 1. 
c 


PHI 

Angle of roll, degrees (<M . 


RE 

Reynolds number (V^D/v). 
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RCORE 


Ratio of the local vortex core radius to 
the local body diameter. Default value is 
.025, and maximum allowable value permitted 
in the code is . 05*SD/r Q . 

XMACH Mach number, M ; 

00 

Item 6 is a single card containing the specification of 
the axial extent of the run and certain parameters associated 
with the vortex wake. Typical values are shown with the 
variables . 


Item Variable 


Description 


6 XI 

XF 
DX 


Initial x-station. XI > 0 
Final x-station. XF >_ XI 

Increment in x for vortex shedding calcula- 
tion. Typical value, DX ~ D/2 when 
FDX = 0.0. For variable DX, set DX = 0.0 
and FDX > 0.0. 


EMKF 


RGAM 


Minimum distance of shed vortex starting 
position from body surface. 

=1.0, vortices positioned such that 

separation point is a stagnation 
point in the crossflow plane 
>1.0, minimum radii away from body surface 
for shed vortices. Typical value, 
EMKF =1.05 

Vortex combination factor. 

=0.0, vortices not combined; preferred 
>0.0, radial distance within which 

vortices are combined. Typical 
value, RGAM = 0.05 D 
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VRF 


Vortex reduction factor to account for 
observed decrease in vortex strength. 

=0.6, for subsonic flow 
=1.0, for closed bodies, or supersonic 
flow? preferred 

FDX Factor for variable x-increments . When 

FDX > 0.0, DX = FDX * (local body radius). 

DX should be input as zero when this option 
is used. 

Item 7 contains only one variable that is of general 
use in program NOZVTX, and that is the integration error 
tolerance, E5. The variable XTABL is for diagnostic 
purposes only. The next three variables concern the use of 
forced asymmetry for bodies at very high angles of attack. 

The variables are defined as follows. 


Item Variable Description 

7 E5 Error tolerance for vortex trajectory 

calculation. (Relative error in vortex 
position.) Typical range, E5 = 0.01 to 
0.05. 


XTABL 


XASYMI 


XASYMF 


x-location after which a table of corres- 
ponding points between the real body and 
the transformed circle is not printed. 

For diagnostic purposes on noncircular 
shapes. Typical value, XTABL = 0.0. 

Initial x-location at which forced 
asymmetry of separation points is used. 
Typical value, XASYMI = 0.0. 

Final x-location at which forced asymmetry 
of separation points is used. Typical 
value, XASYMF - 0.0. 
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DBETA Amount of forced asymmetry for separation 

points on body, degrees. Typical value, 
DBETA = 0.0. 

Items 8, 9, 10, and 11 provide a table of geometry 
characteristics that must be input for all bodies, circular 
or noncircular in cross section. 


Item Variable 


Description 


8 NXR 

9 XR 

10 R 


11 DR 


Number of entries in body geometry table 
(1 <_ NXR £ 50) . 

x-stations for geometry table (NXR values, 
8 per card) . 

r Q , circular body radius at x-stations, or 
radius of tranformed circle (NXR values, 

8 per card) . For an ellipse, 
r 0 - (a + b)/2. 

dr Q /dx, body slope of transformed shape 
at x-stations (NXR values, 8 per card) . 


Items 12 and 13 are included only if NCIR = 1? that is, 
for elliptical bodies. 


Item Variable 


Description 


12 AE 

13 BE 


a, horizontal half-axis of elliptic cross 
section (NXR values, 8 per card) . 

b, vertical half-axis of elliptic cross 
section (NXR values, 8 per card) . 
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Items 14 through 22 are included only if NCIR s =2; that 
is, for a body with arbitrary cross section that must be 
handled with the numerical transformation and conformal 
mapping procedures. 


Item Variable 
14 MNFC 


MXFC 


Description 

Number of coefficients describing trans- 
formation of arbitrary body to a circle 
(1 £ MNFC £ 25) . 

Number of x-stations at which transforma- 
tion coefficients are defined (1 <_ MXFC <_ 10) 
For a similar shape body at all cross 
sections , MXFC = 1 . 


15 XFC x-stations at which transformation coef- 

ficients are defined. For a similar cross 
section body, MXFC=1 , XFC(l) <_ XI (MXFC 
values, 8 per card) . 


Items 16 and 17 are included when NCF = 0. This block 
of cards is repeated MXFC times, once for each x-station 
input in Item 15. 


Item Variable 


Description 


16 NR Number of coordinate pairs defining the 

body cross section at the axial station 
defined by XFC(J) (2 <_ NR <_ 30) . 


17 


XRC (I , J) 
YRC ( I , J ) 


Coordinates of right-hand side of body. 

The convention for ordering the coordinates 
from bottom to top in a counter-clockwise 
fashion, as shown in Figure 10, is observed 
Right/left body symmetry is required in 
the numerical mapping. (NR cards with one 
set of coordinates per card.) 
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Items 18 through 12 are input in lieu of Items 16 and 17 
when NCF > 0. The values of these variables must be obtained 
from a previous run of the numerical mapping portion of the 
code. The only purpose for providing the option to input the 
mapping parameters is to eliminate the need to recompute the 
numerical mapping in subsequent runs and save computation 


time . 
each 

This set 
x-station 

of items is repeated MXFC times, once for 
input in Item 15 . 

Item 

Variable 

Description 

18 

ZZC ( J) 
RZC (J) 

Axis shift and radius of mapped cross 
section (one card with two values) . 

19 

NR 

Number of sets of body coordinates at a 
particular axial station. 

20 

XRC , YRC 

Coordinates of right side of body defining 
cross section starting at 3 = 0° and 
ending at 0 = 180° (NR cards with one set 
of coordinates per card). See Figure 10. 

21 

THC 

Angle of defined body points in circle 
plane (NR values, 8 per card) . 

22 

AFC 

Mapping coefficients (MNFC values, 6 


per card) . 

Items 23 and 24 are included only if NXFV > 0. Item 23 
specifies the axial stations at which additional output or 
plots are requested. Item 24 is included only if NFV > 0. 
Each card of this item contains the nondimensionalized 
y , z-coordinates of a field point at which the velocity field 
is calculated at each axial station specified in Item 23. 
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Item 


Variable 


Description 

23 XFV x-station at which field point velocities 

are calculated or at which additional out- 
put is printed (NXFV values) . 

Omit Item 24 if NFV = 0. 

24 YFV,ZFV y/r Q , z/r Q -coordinates of field points at 

which velocity field is calculated, 
expressed as a fraction of the local body 
radius (NFV cards with one set of coordi- 
nates per card) . It is important that the 
field points lie outside the body surface. 

Items 25, 26, 27, 28, and 29 are included for a restart 
calculation only, and the condition NVP + NVR + NVM + NVA ^ 0 
must be satisfied. It is not necessary to input all four 
types of vorticity in a restart calculation. 

Item 25 is a single card containing the nose force and 
moment coefficients at the restart point; that is, X = XI . 
These values may be set equal to zero, but the resulting 
forces and moments from the current calculation will be in 
error . 


Item Variabl e 
25 CN 

CY 
CM 
CR 
CSL 


Description 

Normal-force coefficient. 
Side-force coefficient. 
Pitching-moment coefficient. 
Yawing-moment coefficient. 
Rolling-moment coefficient. 
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Item 26 is a block of NVP cards to specify the positive 
separation vorticity on the right side of the body. Omit 
this item if NVP = 0. The variable XSHEDP associated with 
each vortex is used only to identify individual vortices 
and permit the user to follow individual vortices during the 
calculation . 


Item 

Variable 

Description 

26 

GAMP 

r/V^, positive separation vorticity on 
right side of body. 


YP,ZP 

coordinates of discrete vortices on right 
side of body at starting point (XI) . 


XSHEDP 

x-location at which individual vortex was 
shed (may be 0.0) . 

(Item 

26 consists 

of NVP cards.) 

Item 27 is a block of NVR cards to specify the second- 
ary or additional vorticity on the right side of the body or 
at any other position in the field. This array is a con- 
venient way to put an arbitrary cloud of additional vorticity 
which is to be maintained separately from the other body 
vorticity in the field. Omit this item if NVR = 0. 

Item 

Variable 

Description 

27 

GAMR 

r/V^, reverse flow or additional vorticity 
on right side of body. 


XR, ZR 

coordinates of discrete vortices on right 
side of body at starting point (XI) . 


XSHEDR 

x-location at which individual vortex was 
shed (may be 0.0). 

(Item 

27 consists 

of NVR cards . ) 
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Item 28 specifies the negative separation vorticity on 
the left side of the body, analogous to Item 26. Omit this 
item if NVM = 0 or if ISYM =0. In the case of a symmetric 
flow field (ISYM=0) , this block of vorticity is automatically 
defined as the mirror image of the positive vorticity input 
in Item 26. 


Item Variable Description 

28 GAMM r/V^, negative separation vorticity on left 

side of body. 


YM,ZM 


Coordinates of discrete vortices on left 
side of body at starting point (XI) . 


XSHEDM x-location at which individual vortex was 

shed (may be 0.0). 


(Item 28 consists of NVM cards.) 


Item 29 specifies the secondary or additional vorticity 
on the left side of the body or at any other position in the 
field. This block of vorticity is analogous to Item 27, and 
it is omitted if NVA = 0 or if ISYM =0. As with the previous 
item, this block of vorticity is automatically defined as 
the mirror image of the negative vorticity input in Item 27 
if symmetry is required by ISYM = 0. 


Item Variable Description 

29 GAMA r/V^, reverse flow or additional vorticity 

on left side of body. 


YA, ZA 


Coordinates of discrete vortices on left 
side of body at starting point (XI) . 


XSHEDA x-location at which individual vortex was 

shed (may be 0.0). 


(Item 29 consists of NVA cards.) 
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The above 29 items make up the body vortex shedding 
portion of the input deck. 

Items 30 through 43 are associated with the paneling 
of bodies with noncircular cross sections and are taken from 
the input specifications for program WDYBDY (Ref. 10) .. 

Item 44 is used to describe bodies with circular cross 
sections. As described above, it was desired to keep the 
panel code input as unchanged as possible from the original 
WDYBDY code; however, the resulting input deck requires the 
specification of several variables over which the user has 
no control. The input deck considered herein contains only 
those variables which are optional to the user. Other 
nonoptional variables have been hard coded with their 
appropriate value, and these should not be changed by the 
user . 

Omit Items 30 through 43 if NCIR = 0. 

Variable Description 

TITLE1 One card of identification information to 

be printed ahead of the paneling output. 

Item 31 is a single card containing two indices, 
each index is right justified in a three column field. The 
purpose of IXZSYM is to define the manner in which the body 
paneling is carried out for different flow conditions. A 
symmetric body in a symmetric flow ( 4> = 0°) can benefit 
from the IXZSYM = 0 option, which panels only one half of the 
body and then utilizes the symmetry characteristics to reduce 
computation time. The same symmetric body in an asymmetric 
flow condition (<f> f 0°) must be paneled in total, IXZSYM = 1, 
to pick up the nonsymmetry in the loading around the body. 

The body symmetry is used to reduce the input required as 
described later. 


Item 

30 
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Item Variable 


Description 


31 IXZSYM Panel symmetry option. 

=0, panel symmetric half of body for 
symmetric flow (<j> - 0°) 

-1, panel complete symmetric body using 
+y side geometry for a symmetric body 
in asymmetric flow (<j> ^ 0°) 

=-l , panel complete configuration using 
geometry of both sides 


ITBSYM -0, configuration has top/bottom symmetry 

=1, no top/bottom symmetry (e.g., see 
Fig. 10) 


Item 32 contains indices to specify the body geometry. 


Item 

Variable 

Description 

32 

J2 

Body geometry specif ication option. 

=0, no body (not recommended) 

=1 , geometry for arbitrary shaped body 
=-l, circular body defined by cross- 


sectional area at XFUS stations (not 
recommended) 

=-2, circular body defined by radius at 
XFUS stations 

=-3, elliptic body defined by both semi- 
axes at XFUS stations 

J6 =0, cambered body - not available 

=1, body symmetrical with respect to 
XY-plane (uncambered body) 

=-1, circular or elliptic body with 
J2 ^ 0 (preferred value) 
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NRADX 


Number of points used to represent the 
body segment about the circumference. If 
the configuration is symmetric 
(IXZSYM = 0, 1) , NRADX is input for the 
half section. If the entire configuration 
is input (IXZSYM = -1) , NRADX is input for 
the full section. The half section or 
full section is divided into NRADX-1 equal 
angles, and the y- and z-coordinates of 
the panel corner points are defined by the 
intersections of these meridian angles 
with the body surface (3 £ NRADX £ 20) . 

It is the experience of the authors 
that the best panel representation of the 
body is obtained when NRADX is an even 
number; that is, there should be an odd 
number of panels on the half-body. It is 
very important for the pressure calcula- 
tion that there be a panel control point 
on the 3 = 90° meridian. 

NFORX Number of axial body stations used to 

define geometry (2 <_ NFORX <_ 3 0) . 

The next block of input specifies the geometry of the 
body for purposes of laying out the panels. This input is 
somewhat redundant with the input in Items 8 through 13 ; 
however, it is important that the paneling geometry be 
separate from the previous geometry but consistent with it. 
In many cases, the paneling geometry may be specified in a 
more coarse grid than that required by the separation 
portion of the program. 
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Item 33 is a block of cards containing the axial stations 
at which body geometry is input for paneling purposes. The 
values are input ten to a card in seven column fields with a 
decimal point required. 

Item Variable Description 

33 XFUS x-stations at which body paneling geometry 

is input (NFORX values) . 

Item 34 is inlcuded only if J2 = -1 (Item 32) . The 
circular body cross-sectional areas are specified at the 
XFUS stations defined in Item 33. 

Item Variable Description 

34 FUSARD Circular body cross-sectional areas at 

XFUS stations (NFORX values) . 

Item 35 is included only if J2 = -2 (Item 32) . The 
circular body radii are specified at the XFUS stations 
defined in Item 33. 

Item Variable Description 

35 FUSRAD R, circular body radii at XFUS stations 

(NFORX values) . 

Items 36 and 37 are included only if J2 = -3 (Item 32) . 
The elliptic body horizontal and vertical semi-axes are 
specified at the XFUS stations defined in Item 33. 
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Item 


Variable 


Description 


36 FUSBY Elliptic body horizontal semi-axes at 

XFUS stations (NFORX values) . 

37 FUSAZ Elliptic body vertical semi-axes at 

XFUS stations (NFORX values) . 

Item 38 is included only if J2 = 1 (Item 32) . The 
y , z-coordinates of points on an arbitrary cross section 
body are specified in Item 38, and there is one set of 
Item 38 for each of the axial stations defined by XFUS in 
Item 33. The convention for ordering the coordinates from 
bottom to top is observed (e.g., Fig. 10). If the full 
cross section is to be specified (IXZSYM = -1, Item 31) the 
ordering continues counter-clockwise from the top of the 
body back to and including the bottom point to close the 
cross section. 

Item Variable Description 

38 YJ , Z J y , z-coordinates of arbitrary body at XFUS 

stations (NRADX cards) . 

Item 38 is repeated NFORX times, once for each XFUS station. 

The following portion of the input deck is described as 
Auxiliary Input Cards in Appendix K of Reference 10. The 
purpose of this input is to provide a simple means to adjust 
the actual panel layout keeping the geometry specified above 
unchanged . 

Item Variable Description 

39 TITLE2 One card of identification information for 

paneling geometry. 
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Item 40 contains an index to specify the output for the 
body panel portion of the program. Because of the quantity 
or output possible, it is recommended that |lPRINT| £ 1 be 
the highest level of output specified. Note that the 
values of the first two indices on this card are fixed. 

Variable Description 

IPRINT Panel output option. 

-0, minimum output (recommended) 

=1 , output panel corner coordinates, final 
source strengths, and velocities at 
panel control points (recommended) 

- 2 , same as IPRINT=1 and source strengths 
from previous iteration (not 
recommended) 

=3 , print source strength for each itera- 
tion (not recommended) 

=4 , same as IPRINT=3 and complete aero- 
dynamic influence matrix (not 
recommended) . 

If IPRINT < 0, the panel geometry is included with the output 
specified by | IPRINT | . This option provides a useful means 
to check the paneling arrangement. 

Item 41 is a single card containing two indices in three 
column fields. These indices provide a conventional means to 
change the panel layout without changing the geometry input. 

Item Variable Description 

41 KRADX Number of meridian lines used to define 

panel edges on the body segment (see note 
for NRADX in Item 32) . 


Item 

40 
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KRADX 


=0, meridian lines are defined by NRADX 
in Item 32 

>0, meridian lines are calculated at 
KRADX equally spaced points 
(3 £ KRADX £ 20) 

<0, meridian lines are calculated at 

specified values of PHIK in Item 42 
(3 < | KRADX | £ 20) 

For symmetric configurations (IXZSYM= 0,1) , 
KRADX is the number of meridians on the 
half section. For full configurations 
( IXZSYM = -1) , KRADX is the number of 
merdians on the full section including 
meridians at 0° and 360°. 

KFORX Number of axial stations used to define 

leading and trailing edges of panels on 
the body . 

=0, the panel edges are defined by NFORX 
and XFUS in Items 32 and 33, 
respectively 

>0, the panel edges are defined by XFUSK 
in Item 4 3 (2 <_ KFORX £ 3 0) 

Item 4 2 is included if KRADX <_ -3 in Item 41. This 
item provides a means to easily modify the panels by changing 
the body merdian angles which specify the panel corner 
points. The purpose of this option is to allow the user to 
pack additional panels into a region of the body which is 
changing rapidly. 

Item Variable 
42 PHIK 

51 


Description 

Body merdian angles (degrees) with 
PHIK = 0° at the bottom of the body and 
PHIK = 180° at the top ( | KRADX | values). 



Item 43 is included only if KFORX is nonzero. The 
variable XFUSK defines the axial positions of the panel edges. 
If this item is not included, the panels are specified by 
XFUS in Item 33. The purpose of this item is to simplify 
the changing of a panel layout without changing the entire 
geometry input deck. 

Item Variable Description 

43 XFUSK Axial stations of leading and trailing 

edges of body panels (KFORX values) . 

Item 44 is included only if NCIR = 0; that is, for 
bodies with circular cross sections. 


Item Variable Description 

44 NXBODY Number of line sources/sinks and line 

doublet singularities distributed along 
body centerline (1 £ NXBODY £ 100) . 
Default value is NXBODY = 51. 


NCODE 


XNOSE 


XLBODY 


Integer control index for specifying fore- 
body shape over length of XNOSE. 

NCODE =0, Parabolic 

=1, Sears-Haack 
=2, Tangent ogive 
=3, Ellipsoidal 
=4 , Conical 

Length of nose part of body, measured from 
nose tip. 

Length of body. 
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This concludes the input deck for a single run of program 
NOZVTX. The program is not set up to stack input for 
multiple cases because of the sometimes long and generally 
unpredictable run times. 

Input Preparation 

In this section, some additional information is provided 
to assist the user in the preparation of input for certain 
selected problems. The previous section on the input 
description must be used to understand the individual vari- 
ables which go into NOZVTX, and this section will permit the 
user to select the appropriate input to get optimum results 
from the code. 

Panel layouts .- The choice of panel arrangement to 
represent the noncircular body shape is perhaps the most 
important part of the input selection as the quality of the 
panel layout can determine the quality of the predicted 
results. Simply using the largest possible number of panels 
is one solution; however, this will not guarantee an 
optimum panel layout, nor will it guarantee the best results. 
In the following paragraphs, some guidelines are presented 
for various situations to optimize both results and computa- 
tion time whenever possible. 

The elliptic cross section will be considered first 
because the code has been checked out for a number of 
different elliptic bodies. In each case to follow, right/ 
left symmetry of the panels (IXZSYM = 0 or 1 in Item 31) is 
assumed; therefore, only one side of the body will be shown 
to illustrate the panel layout. 

A 2:1 elliptic cross section is shown in Figure 11 where 
several different panel arrangements have been selected. In 
Figure 11(a), nine program-determined panels (KRADX = 10 
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in Item 41), each specified by a 20° central angle, are shown. 
The only problem area is near the shoulder of the ellipse 
where the outer three panels miss the actual surface by a 
considerable amount. Increasing the number of panels to 
10 (KRADX = 11) produces the arrangement shown in 
Figure 11 (b) . Even though an additional panel has been 
included, the region near the shoulder is not modeled 
better than before. In fact, this layout is missing an 
important physical effect because no panel has the same 
orientation as the shoulder (dz/dy = °°) , and the entire 
effect of the shoulder is lost. 

An improved panel layout is shown in Figure 11 (c) where 
nine panels are specified by their meridian angles 
(KRADX = -10) . In this layout, the panels are concentrated 
in the region of highest curvature near the shoulder and 
spread out over the lessor curved region. Note that an 
odd number of panels is specified so that the outer panel 
nearest the shoulder has the proper slope. Top and bottom 
symmetry of the panels is shown, and this is a recommended 
arrangement for cross sections with top/bottom symmetry. 

For this particular cross section, additional panels can 
be specified for better resolution at the cost of increased 
execution time. It is important that the user check the 
panel layout used by the code to ensure reasonable modeling 
of the actual surface. 

Different panel arrangements on a 3:1 elliptic cross 
section body are shown in Figure 12. The sharper curvature 
of the shoulder of this shape creates a bigger modeling 
problem than the 2:1 ellipse shown previously. In 
Figure 12(a), a 13-panel arrangement (KRADX = 14) with equal 
spacing is shown, and it is obvious that a large portion 
of the shoulder region is not modeled properly for the 
prediction of flow field characteristics near to the body. 
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For this particular shape, more panels are required for 
better resolution. Additional panels are added in the 
region 69° < 3 < 111° to produce the 19-panel layout 
(KRADX = -20) shown in Figure 12(b) . 

An example of reasonable panel arrangements on a lobed 
cross section body without top/bottom symmetry is shown 
in Figure 13. The shape shown in this figure is that 
described in References 17 and 18. The panel edges do not 
lie precisely on the actual body shape, a problem which is 
discussed in the numerical mapping section to follow. The 
9-panel layout in Figure 13 (a) is a resonable representation 
for preliminary calculations; however, certain detail is 
missing in the regions of greatest curvature. Additional 
panels are included in these regions, and the resulting layout 
is shown in Figure 13 (b) . Custom panel layouts are shown 
for both lobed-body examples to take advantage of the fewer 
panels required on the flat portion of body. 

It is difficult to specify rigid guidelines for panel 
layout on arbitrary cross sections. The examples shown 
above were determined by this author after several trials, 
and proper panel selection becomes easier as the user gains 
experience with the code. The user should remember that 
linear interpolation of panel characteristics is used between 
panel control points in the current version of the code. 

Based on the experience of the authors, it is always useful 
to compare the paneling layout with the actual body shape, 
and generally, good engineering judgement will allow the 
user to select a reasonable panel layout. 

Nothing has been said about axial panel layout (NFORX 
or KFORX) in the above discussion. Axial spacing can be 
critical to the predicted results, and as mentioned above, 
good engineering judgement will usually produce an adequate 
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layout. Linear interpolation is also used in the axial 
direction to determine panel characteristics between panel 
control points. This procedure can be used to great 
advantage for a specific class of body, the cone. Since the 
predicted linear results on a cone are the same at each 
axial cross section, fewer panels are required in the axial 
direction. The major requirements are the following. 

1. There should be one row of panels sufficiently 
close to the nose tip or starting point (XI) for the 
calculation that a control point can fall between the 
starting point and the nose tip. 

2. There should be a row of panels with control points 
located just aft of the last axial station of the calcula- 
tion (XF) . 

3. If a conical nose is attached to the cylindrical 
afterbody, there should be enough axial stations specified 
near the junction of the nose and afterbody that adequate 
body characteristics are available for interpolation 
purposes . 

An example of a panel arrangement for a cone forebody 
is presented in the sample cases. 

Numerical mapping .- The specification of the appropriate 
numerical mapping parameters (Items 14 through 22) depends on 
the shape of the cross section. In the interest of optimum 
computation time, the fewest possible terms in the transforma- 
tion series (MNFC in Item 14) should be used which will 
produce an adequate mapping. Each different shape should be 
checked by the user to determine the proper number of 
coefficients. For example, a simple cross section which is 
similar to a circle may require as few as 10 coefficients. 

The shape shown in Figure 13 required 20 coefficients for 
an excellent mapping. If a large number of calculations are 
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to be made for a noncircular shape, it is worthwhile for the 
user to try different numbers of coefficients in an effort 
to find an optimum number. Not only should the shapes be 
compared carefully, the unseparated flow pressure coefficients 
should be compared to evaluate the quality of the mapping. 

The numerical mapping now included in NOZVTX has been 
used for a variety of cross section shapes, and very few 
problems have been observed. One general problem area is 
bodies with concave sides. The lobed body in Figure 13 
has a moderate concave section on each side, and the mapping 
procedure had difficulty converging on a set of coefficients. 
When the body shape was modified slightly, as shown in 
Figure 13(a) , to give the concave region a flat or slightly 
convex shape, the mapping converged easily. This small 
change in the shape probably has very little effect on the 
predicted pressure distribution. This is an example of the 
type of problems to which the user should be alert. 

Integration interval .- Choice of the appropriate 
integration interval (DX in Item 6) depends to a large extent 
on the type of results of interest. The code permits the 
interval to be a constant (FDX = 0 in Item 6) or a function 
of the local radius (FDX > 0) . If the user is interested 
in detailed nose loads and separation on the nose, where 
the radius is changing rapidly, the latter option should be 
used with 0.5 £ FDX <_ 1.0. This option can lead to a wide 
range of vortex strengths which, on occasion, has been 
determined to be the cause of certain vortex tracking 
difficulties. The user should be alert to this possibility. 
The constant interval option (DX = constant) is the best 
choice when the results of a major interest are the vortex 
clouds near the end of the body. 
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Vortex core .- The core size specification is referenced' 
to the local body diameter (RCORE in Item 5) . The default 
value is 0.025, and this value is adequate for most cases. 

The larger the core radius specified, the more the vortex 
effects are reduced. To eliminate the possibility of 
negating the vortex effect entirely, a maximum vortex core 
of .05 times the maximum body diameter is built into the 
code. Except in unusual situations, the default value is 
recommended, and this is the value used in the comparisons 
with data shown in a later section. 

Sample Cases 

In this section, 11 sample input decks are described 
to illustrate the input preparation and use of the program 
for a wide range of configurations and flow conditions. 
Insofar as possible, each major option in the program is 
illustrated in the set of sample cases. The major features 
of each sample case are identified in Table II to assist the 
user in finding an appropriate example to suit a specific 
need . 

Sample Cases 1 through 11 shown in Table II are designed 
to exercise the most common options of program NOZVTX, and 
Sample Case 1 is included specifically to provide a short 
computer run for testing purposes. The input decks for the 
sample cases are shown in Figures 14 (a) through (k) . The 
sample cases are based on typical bodies used in the data 
comparisons to follow in a later section, and the reference 
in which the configuration is defined is noted on each 
sample case. 

Sample Case 1 is a test case based on the 3:1 elliptic 
missile body in Reference 19; however, only a portion of 
the nose is considered to keep the execution time and output 
at a minimum. The purpose of this run is to provide the 
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user with a short run to check out the code. The output from 
this case is described in the next section. Sample Case 2 
in Figure 14(b) is a restart of the previous run, and it 
illustrates the continuation of a calculation. Sample 
Case 3 in Figure 14 (c) is the input deck for a complete run 
of this 3:1 elliptic missile. The first three sample cases 
do not necessarily use the optimum panel arrangement in an 
effort to minimize execution time. For production runs, the 
authors recommend a panel layout similar to that shown 
in Figure 12 (b) . 

Sample Case 4, Figure 14(d), illustrates the missile 
body in a rolled condition (<J> = 45°) . The lack of flow 
symmetry and consequently the lack of right/left panel 
symmetry dictate the use of a smaller number of panels 
because the entire cross section must be paneled. This 
run is typical of any noncircular body in a rolled flow 
condition . 

Sample Cases 5 and 6, Figures 14(e) and (f) , illustrate 
input decks for circular cross section bodies (Refs. 20 
and 21) . Since axis singularities are used in place of 
panels, the input is shorter and easier to prepare. The 
primary objective of Sample Case 6, an axisymmetric body 
at a c = 0° , is the prediction of the surface pressure 
distribution with no separation effects included. 

Sample Cases 7 and 8, Figures 14(g) and (h) , are input 
decks for a 2:1 elliptic cross section forebody from 
Reference 18. The first case is for a standard panel 
layout determined by the code. The panels for the second 
case are tailored to the cross section shape by specifica- 
tion of the central angles of the panel edges. 

Sample Case 9 is a 3:1 elliptic cross section conical 
forebody (Ref. 17) . The input specifies a variable axial 
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integration interval, and the axial panel distribution is 
determined by the conical shape. Because of the character 
of conical flow, a very small number of axial panels is 
required as discussed in a previous section on input 
preparation. 

Sample Case 10, Figure 14(j), is a lobed cross section 
forebody from Reference 18. This is an arbitrary noncircular 
body and numerical mapping is required for solution as 
indicated in the input. Since the cross section is changing 
at each x-station, the cross section shape must be input at 
a number of stations thus accounting for the large quantity 
of input. The number of axial stations must be sufficient 
to provide acceptable accuracy in the body shape using linear 
interpolation between input stations. 

Sample Case 11, Figure 14 (k) , is set up for the 2:1 
elliptic cross section missile body in Reference 22. The 
panel layout for this body is the same as that shown in 
Figure 11(c) . 


Output 

A typical set of output from program NOZVTX is described 
in this section. In general, the output quantities are 
labeled and each page is headed with appropriate descriptive 
information. The actual output obtained is determined by the 
output options selected by the user. For purposes of this 
description, representative pages from Sample Case 1, 
described in a previous section, are presented in Figure 15. 

The total output from Sample Case 1 is contained on 30 
pages, including pages which are continuations of another 
page. The first page shown in Figure 15(a) is headed by 
the title cards input as Item 3. The remainder of this 
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page, the following two pages shown in Figures 15(b) and (c) , 
and the top of the fourth page shown in Figure 15 (d) contain 
all the input data with appropriate labels. Body source 
panel geometry is shown on pages 5 through 13, Figures 15(e) 
through (m) . Page 5, shown in Figure 15(e) and the following 
two pages shown in Figures 15(f) and (g) indicate the panel 
corner coordinates. Pages 8, 9, and 10, shown in 
Figures 15(h) , _ (i) , and (j) , contain a table of the body 
panel centroid coordinates. This is followed by a table of 
body panel areas and inclination angles. The inclination 
angles DELTA and THETA orient the panels in the body 
coordinate system. Detailed descriptions of these angles are 
given in Reference 10. This table continues onto page 12 
and the top of page 13, shown in Figures 15 (£) and (m) . The 
next item on page 13, GB(N), is the list of source strengths 
for the body at angle of attack. This is followed by a table 
of velocity components at the control point of each panel. 
Each component is labeled, and the table continues onto pages 
14 and 15 shown in Figures 15 (n) and (o) . The last item on 
page 15, under the heading MODIFIED SOURCE STRENGTH, is 
the source strength on each panel of the body with no cross- 
flow influence. This panel strength represents the volume 
solution of the body for purposes of tracking vortices. 

This list of source strengths continues onto page 16 shown 
in Figure 15 (p) . 

The next two pages, shown in Figures 15 (q) and (r) , 
contain the predicted pressure distribution around the body 
at the first axial station. A summary of the geometric 
characteristics and the Mach number is printed at the top 
of page 17, followed by the pressure distribution. The 
columns are labeled, but a brief explanation of the printed 
items follow. The first three columns, Y, Z, and BETA, 
locate the specific points on the body cross section. U/VO 
is the local axial velocity referenced to the free-stream 
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velocity. V/VO and W/VO are the local lateral and vertical 
velocities, respectively. These velocities include all 
perturbation effects including vortex induced velocity com- 
ponents. The column VT/VO is the total velocity in the 
crossflow plane given as the vector sum of the V/VO and W/VO 
velocities. A positive VT/VO represents flow in a counter- 
clockwise direction when viewed from the rear. CP (M) is the 
compressible pressure coefficient as calculated from 
Equation (29) . DPHI/DT is the unsteady part of the 
incompressible pressure relation as represented by the last 
term in Equation (31) . CPZ is the steady part of 
Equation (31), and CP (I) is given by Equation (31) in its 
entirety. VCP/VO is the local two-dimensional crossflow 
plane velocity whose positive sense is the same as that for 
VT/VO . 

The local force and moment coefficients are printed at 
the middle of page 18 . CN (X) and CY (X) are local force 
coefficients calculated by Equations (39) and (43) , 
respectively. The remainder of the forces and moments 
represent total quantities on the nose between the origin 
and the local axial station. The normal force and pitching 
moment coefficients, CN and CM, are calculated with 
Equations (40) and (41), and the axial center of pressure 
of the normal force, XCPN , is given by Equation (42) . XCPN 
is not normalized by the reference length on page 7. The 
corresponding side force and yawing-moment coefficients and 
center of pressure of the side force are given by 
Equations (44) , (45) , and (46) , respectively. The rolling- 
moment coefficient, CSL, is positive in a counterclockwise 
sense when viewed from the rear of the body. The small 
value shown in this output is due to numerical round-off 
error during the calculation. The same is true for all 
lateral coefficients which appear as very small, nonzero, 
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values when a symmetric body and symmetric flow conditions 
dictate that they are identically zero. 

Page 19 of the output from Sample Case 1, Figure 15 (s), 
summarizes the important points in the pressure distribution 
and the separation calculation. The last two lines of output 
contain specific information on the vortices shed at this 
station. The vortex numbers, NV, assigned here are the 
cumulative numbers of vortices in the flow field. The 
vortex strength, GAM/V, is obtained from Equation (38) where 
the boundary layer edge velocity u /V to is shown as VT/V in 
this line of output. M(K) is the actual distance of the 
vortex from the body surface. The position of the vortex 
in the body plane is given by Y , Z, and BETA. The position 
of the vortex in the circle plane is given by the coordinates 
YC , ZC. The radial position of the vortex in the transformed 
circle plane is RG/R. The latter quantity is useful in 
providing a quick estimate of the relative positions of new. 
vortices being shed into the wake. 

Page 20, shown in Figure 15 (t) , summarizes the vortex 
field after the trajectories have been calculated over one 
axial interval. The x-station at the top of this page is 
exactly DX downstream of the previous x-station. Information 
to locate each individual vortex in both the body plane and 
the circle plane is provided. One additional piece of 
information for each vortex is XSHED , the x-station at which 
the individual vortex originated. This is printed to allow 
the user to keep track of any individual vortex throughout 
the trajectory calculation. Each vortex maintains its 
individual identity and the right and left side vortices are 
grouped together on this page. Centroids of the individual 
groups of vorticity also show on page 20. The only exception 
to the above discussion of individual vortices occurs when 
vortices are allowed to combine through use of the parameter 
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RGAM in Item 6 of the input list. When vortices are combined, 
the particular vortices are identified in a special line of 
output . 

This completes one cycle of output. The next page of 
output is the pressure distribution at the new axial station, 
and the format is the same as page 17 shown in Figure 15 (q) . 

At each axial station, the output repeats through the cycle 
shown in Figures 15 (q) , (r) , (s) , and (t) . At specified 

axial stations, certain additional output is available. A 
simple plot of the body cross section and the vortex field 
at x = 4.2 is shown in Figure 15 (u) . This particular output 
is on page 23 of the actual output from Sample Case 1. Each 
type of vortex is shown with a different symbol in this figure. 
Posti've body vorticity shed from the right side of the body 
is shown as an X, and negative vorticity from the left side 
is represented by an 0. If it is requested, the calculated 
velocity field appears as shown in Figure 15 (u). This 
output appears on page 24 of the actual output. The 
coordinates and velocities shown on this page are in the 
body plane. The normalizing diameter, D, in the last two 
columns is the equivalent diameter of the noncircular body. 

The last page of output is shown in Figure 15 (w) . The 
top portion of the output is the summary page for the separa- 
tion calculation at the last x-station, X = XF . This part 
of the output is similar to that shown in Figure 15 (s). 

The lower part of the page is printed only at the completion 
of the calculation, and it contains a summary of the forces 
and moments and the vortex strengths and centroids at the 
end of the body. 
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RESULTS 


For purposes of evaluating the accuracy and range of 
applicability of the engineering prediction method and 
associated program code NOZVTX described above, comparisons 
of measured and predicted aerodynamic characteristics of 
various configurations have been made. The objective of 
the method is to predict the characteristics for noncircular 
cross section bodies, but selected data for bodies of revolu- 
tion will also be considered to illustrate certain features 
of the prediction method. Typical results from the predic- 
tion method are presented in the following sections. 

Circular Bodies 

The prediction method was applied to an ogive-cylinder 
model (Ref. 21) in supersonic flow to evaluate the ability 
of the prediction method to calculate pressure distributions. 
The configuration has a circular cross section and a three 
caliber ogive nose followed by 3.67-caliber cylindrical 
afterbody. Circumferential pressure distributions at a large 
number of axial stations are available for a range of angles 
of attack and supersonic Mach numbers. The free-stream 
Reynolds number, based on diameter, is 0.5 x 10^ for the 
results to follow. 

The model is represented by 51 line sources on the 
body body axis for the calculation at a c = 0° . The measured 
and predicted pressure distributions at three Mach numbers 
are compared in Figure 16 where the multiple data symbols 
at each axial station represent the range of scatter in the 
data. The purpose of this comparison is to verify the body 
model and the pressure calculation technique under conditions 
in which there are neither crossflow nor separation effects. 
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In Figure 16(a) , the agreement is very good at M^ = 
except in the region very near the nose. As seen in 
Figures 16 (b) and (c) , the agreement between measured and 
predicted pressures on the nose gets progressively poorer as 
the Mach number increases. It is possible that nonlinear 
effects caused by the close proximity of the shock wave are 
responsible for this disagreement. For purposes of this 
investigation, a small error in the pressure very near the 
nose will have very little effect on the vortex shedding 
from the remainder of the body. As a further check on the 
body model, a panel distribution of 13 rings of 14 equally 
spaced circumferential panels was used for the same calcula- 
tion. The pressure results are virtually identical to those 
obtained from the line singularity model. 

The prediction method was next applied to the same model 
at a = 20° angle of attack and M =1.6. At this angle of 
incidence, the separation vortex wake on the lee side of the 
body is well developed. Measured and predicted pressure 
distributions at two axial stations on the nose and one 
station on the cylinder are compared in Figure 17 . The 
agreement between experiment and theory at x/D ~ 0.8 is very 
good except on the windward side of the body upstream of 
the minimum pressure point. Similar results are illustrated 
in the lower portion of Figure 17 (a) for a station near the 
end of the nose (x/D = 2.8) . The predicted minimum pressure 
point is upstream of the measured value, but the general 
trend of the predicted pressure distribution is correct. The 
slightly irregular behavior of the predicted curve near 
3 = 130° is caused by local interference of the lee side 
separation vortices. 

Measured and predicted pressures on the cylindrical 
portion of the body (x/D = 5.1) are shown in Figure 17(b) . 

A large vortex wake, also shown in the figure, has developed 
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on the lee side, and it has a dominant effect on the pressure 
distribution. The solid curve, representing the predicted 
pressure distribution in the presence of the vortex wake, is 
in good agreement with the measured results near the mini- 
mum pressure point and on the lee side. The predicted 
pressures are lower than those measured on the windward side 
of the body. There is also some disagreement between 
experiment and theory near 3 = 130° where it appears that 
some of the vortex induced effects are missing from the 
predictions. Contrast the solid curve with the dashed curve 
representing a predicted pressure distribution with no 
separation effects included. 

The prediction method was next applied to an ogive- 
cylinder model for which detailed lee side flow field 
measurements are available (Ref. 20) . This model is an 
ogive-cylinder with a two-caliber nose and a thirteen-caliber 
cylindrical afterbody. The tests were conducted at Mach 
numbers of 2 and 3 at a free-stream Reynolds number, based 
on diameter, of approximately 2x 10 . The data consist of 
lee side velocity components and separation vortex strengths 
and positions. Neither pressure nor force and moment data 
are available from these tests. 

The model is represented by 51 line source and doublet 
singularities on the body axis. The predicted vortex cloud 
pattern at x/D = 10 is shown in Figure 18 for M to = 2 and 
a c = 15° . The vortex shedding is symmetric on both sides of 
the body. Also shown in this figure are the separation 
points located at approximately 86 degrees measured from 
the windward stagnation point. The measured and predicted 
vortex strength and centroid of vorticity are shown to be in 
reasonable agreement in this figure. The predicted vortex 
wake is approximately 9 percent stronger than the measured 
value, and the centroid of the predicted vortex cloud is 
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slightly lower and more outboard than the measured results. 
Again, as a further check on the body model, the body was 
represented by 15 rings of 14 equally-spaced circumferential 
panels. The predicted results were nearly identical in each 
case? however, the execution time for the line singularity 
model was only 20 percent of that required for the panel 
solution. 

Measured and predicted body vortex strengths on the 
same ogive cylinder at Mach numbers 2 and 3 and angles of 
attack 10, 15, and 20 degrees are compared in Figure 19. 

At = 2 , the predicted results are slightly higher than 
those measured. At = 3 , the predicted results are lower 
than those measured on the forward portion of the body and in 
good agreement on the rear portion. The slopes of the 
predicted curves are in good agreement with the experimental 
results, an indication that the predicted shedding rate and 
strengths of the shed vortices are correct. The difference 
in magnitude between the measured and predicted vortex 
strength in Figure 19 is attributable in part to the location 
of the onset of separation. If separation is predicted to 
start further aft on the body than the actual start of 
separation, the predicted vortex strength will be less than 
measured. The maximum error in vortex strength exhibited 
in Figure 19 represents an error in the predicted axial 
position of onset of separation of less than one body diameter 
except at the highest angles of attack. 

Measured and predicted downwash velocity fields in the 
plane of symmetry of the model at = 2 and = 15° are 
shown in Figure 20 for three axial stations. The results 
are generally in good agreement, and the predicted velocity 
distribution exhibits the correct trend in every case. 

The vortex field used to calculate the induced velocities 
at x/D = 10 is the distribution shown in Figure 18. 
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Additional comparisons of measured and predicted down- 
wash velocities on the lee side of the body at x/D = 10 are 
shown in Figure 21. In Figure 21(a), the velocities are 
calculated along a line at z/D = .75 which passes through the 
predicted centroid of vorticity. The agreement is only fair 
in the vicinity of the centroid (y/D ~ 0.4), an indication 
that the present discrete vortex cloud model is not completely 
reliable inside the cloud. In Figure 21(b), the measured 
and predicted velocities are compared along a line at 
z/D = 1.05 which passes above both the measured and predicted 
centroid of vorticity. These results indicate the correct 
trend in the predicted velcoity distribution? however, since 
the strength of the wake vorticity is nearly correct, it 
appears that a cloud model that is higher off the body is 
required to bring the predicted velocities into better agree- 
ment with measurements. The need for a higher cloud is 
verified in Figure 18 where it is shown that the actual 
centroid of vorticity is approximately .15D above the pre- 
dicted centroid. As shown in Figure 19, the strength of the 
vortex cloud is in good agreement with experiment. 

Similar comparisons of measured and predicted downwash 
velocity fields on the lee side of the ogive cylinder at 
Moo = 3 , a c = 15° are shown in Figures 22 and 23. The results 
on the plane of symmetry on the model lee side are shown in 
Figure 22 where it is obvious that the comparisons are not as 
good as those shown in Figure 20 for = 2. The general 
trends are correct, but the magnitude of the velocity 
components is in error. This is likely due to a combination 
of circumstances. For example, in Figure 19, the predicted 
vortex strength corresponding to M w = 3, a c = 15° is less 
than that measured at x/D = 7 and in good agreement with 
measurements at x/D » 13. As seen in Figure 22, the vortex 
induced velocity is less than that measured at x/D = 7 and 
nearly correct at x/D = 13. The fact that the predicted 
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velocity profiles are not in agreement may be due, in part, 
to the problem discussed previously, the error in the 
predicted centroid of vorticity. The velocity profiles 
indicate that the predicted centroid of the separation 
vorticity may be too low, a possible problem at higher Mach 
numbers . 

Analogous to the M m = 2 results in Figure 21, comparisons 
of measured and predicted velocity profiles on the lee side 
of the body at x/D = 10 are shown in Figure 23 for = 3. 

The results look very good in Figure 23(a) at z/D = .75 where 
the region of comparison is in the vortex cloud and near the 
vertical centroid of vorticity. In Figure 23(b) , the agree- 
ment between the velocity profiles is poor with the predicted 
velocities being much less than those measured. 

Elliptic Bodies 

The prediction method has been applied to a number of 
different elliptic cross section bodies. A variety of 
results are presented for several configurations under a 
range of flow conditions, the major objective being to 
examine a large number of configurations to identify the 
strengths and weaknesses of the code . 

The 3:1 elliptic cross section missile of primary 
interest in these comparisons is the sharp-nose body described 
in References 19, 23, and 24. The tests were conducted at 

= 2.5 and a free-stream Reynolds number of 2 x 10^ per foot. 
The sharp-nose body considered in the comparisons to follow 
was modeled with 18 circumferential panels around the cross 
section and 29 axial rings of panels. Sample Cases 3 and 4 
described in Table II and Figures 14(c) and (d) , respectively, 
produce the results to follow. 
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Flow visualization via vapor-screen photographs 
(Ref. 23) and pressure data (Ref. 24) are available for 
comparison with the theoretical results. For better 
visualization of the vortex development along the body, the 
photographs were digitized. These results are presented 
in the upper left corner of Figure 24 for the missile at 
= 2.5, - 20°, and <p = 0° and 45°. Cross sections 

through the body at three axial stations are shown. The 
extent of the vortices from the vapor screen photographs 
are shown by the solid curve, and the predicted discrete 
vortex wake is shown by the individual vortices connected 
by a dashed line. The unrolled flow condition (<J> = 0°) 
shown in Figure 24 (a) illustrates very good agreement between 
the predicted vortex cloud and the vapor screen results. 

The strength of the predicted vortex cloud is shown at each 
axial station, and because of symmetry, the left and right 
vortices are of equal magnitude and opposite sign. 

The rolled flow condition (<J> = 45°) is shown in 
Figure 24 (b) . To be consistent with the output from the 
code, the body is shown unrolled with the flow angle rolled 
45 degrees. The general character of the predicted vortices 
is in fair agreement with the experimental results; however, 
there are certain areas of disagreement. The theoretical 
vortex on the right side has rolled up much tighter than 
the vapor screen indicates, and the left theoretical vortex 
has not moved out away from the body as far as the actual 
vortex, nor has it rolled up as tightly as the actual 
vortex . 

Comparison of measured and predicted pressure distri- 
butions at x/L = 0.60 on the body are shown in Figure 25 for 
the unrolled and rolled condition. The unrolled results 
shown in Figure 25 (a) illustrate some of the same problem 
areas discussed for circular bodies. The predicted 
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stagnation pressure on the windward side is in good agree- 
ment with experiment, but the predicted pressure around the 
shoulder region of the cross section is much lower than the 
measurements. The pressure on the lee side is also under- 
predicted, possibly due to vortex interference effects. As 
shown in Figure 24 (a) , there is a strong concentrated vortex 
field developed on the lee side at this axial station 
(x/L = 0.60) . The irregularities in the predicted pressures 
near the shoulder region are caused by local vortex inter- 
ference effects. As noted in the discussion for the circular 
bodies, the discrepancy in the pressure comparisons on the 
windward side of the body near the shoulders may be caused 
by nonlinear effects induced by the close proximity of the 
shock wave. 

The measured and predicted pressures on the 3:1 elliptic 
body in a rolled flow condition are compared in Figure 25 (b) . 
The qualitative agreement is very good, but a problem on the 
windward side of the body is in evidence. The vortex induced 
effects cause the irregularities in the predicted pressures, 
and it appears that the induced effects are too large on the 
positive lee side (0 ~ 320°). In Figure 24(b) it is obvious 
that the close proximity of the predicted wake will have a 
large effect on the pressure distributions, and it is likely 
that the predicted wake is causing the lack of agreement in 
this region. 

NOZVTX was also applied to a conical forebody with a 
3:1 elliptic cross section shape (Ref. 17). The major axis, 
in planform, is defined by a 20-degree half angle, and the 
total length of the forebody is 14 inches. The test was 

g 

conducted at a free-stream Reynolds number of 2 x 10 per foot. 
The panel model representation consists of 17 circumferential 
panels on the half body (34 panels on the total cross section) 
and four axial rings of panels. The limited number of axial 
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panels is permissible because of the similar characteristics 
on a conical configuration in supersonic flow described in 
the previous section on input preparation. The detailed 
panel layout is described in Sample Case 9. 

Comparisons of measured and predicted pressure distri- 
butions on the conical forebody are shown in Figure 26 for 
two axial stations. Except for small differences on the lee 
side of the body, the measured pressure distributions are 
nearly identical at the two axial stations. In Figure 26(a) , 
the theoretical pressure distribution is not in good agree- 
ment with experiment on the windward side of the body, 
particularly near the shoulders. On the lee side, two pre- 
diced curves are shown; the dashed curve represents no 
separation vortex effects, the solid curve has vortex effects 
included. There is considerable difference in these results, 
all of the difference attributable to the vortex induced 
effects. The vortex effects improve the predicted pressures 
in the region around 3 = 120° on the body, but they detract 
from the agreement near the lee side stagnation point, 

150° < 3 < 180°. This may be an indication that the actual 
lee side vorticity is rolled up tightly around 3 — 120°, 
whereas the predicted vortex cloud is distributed over the 
entire lee side of the body as illustrated in Figure 27. 

The prediction method has been applied to the elliptic 
cross section body described in Reference 22 and illustrated 
in Figure 28. The major axis of the cross section is 
specified by an ogive distribution on the nose, and the 
cross section is a 2:1 ellipse over the entire body length. 
The equivalent base diameter of this body is 1.125 inches. 
Circumferential pressure measurements at numerous axial 
stations are available at = 1.5 and a range of angles of 
incidence and angles of roll. The tests were conducted at a 
free-stream Reynolds number, based on the equivalent 

g 

diameter, of approximately 0.3 x 10 . 
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This 2:1 elliptic cross section body is modeled with 
18 circumferential panels around the complete cross section 
and 15 axial rings of panels. No smoothing is used in any 
of the calculations to follow. The general layout of the 
geometry is illustrated in Sample Case 11 in Figure 14 (k) . 

Measured and predicted pressure distributions on the 
elliptic cross section body at M^ - 1.5, a = 15°, and 
several roll angles are compared in Figure 29 at two axial 
stations on the nose . The agreement between experiment 
and theory for cj) = 0° is better away from the tip of the 
nose as shown in Figure 29(a). This is consistent with the 
circular body results described previously. The sparse 
nature of the measurements make it difficult to fully 
evaluate the prediction method other than to observe that 
the basic trends of the pressure distribution are predicted. 
At x/D ~ 2.1, the separation vorticity has grown in strength 
and scope, and the major region of vortex wake influence 
covers the region between & = 90° and 130°. 

The predicted pressure distributions on the same body 
at roll angles of 22.5° and 45° are shown in Figures 29(b) 
and 29(c), respectively. The predicted results are generally 
in good agreement with experiment, and the correct trends of 
the data are modeled by the theory in each case. The local 
irregularities in the predicted curves are due to the 
influence of the separation vortices. 

NOZVTX was next applied to a 2:1 elliptic cross section 
cone (forebody 4 in Ref. 17) in an unrolled flow condition. 
The pressure results for - 1.7 and a c = 20.35° are shown 
in Figure 30 where it is obvious that agreement is poor 
over much of the surface. The lack of agreement on the wind- 
ward side and near the body shoulder is typical of previous 
results on elliptical cross section bodies. On the lee side, 
there is a large effect of the predicted vortex field which 
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is spread out over much of the lee side. The effects of the 
vortices on the measured pressures appear to be more limited 
in scope. These results are the same as those shown on the 
3:1 elliptic cone forebody in Figure 2*6 (a) . 

The final 2:1 elliptic cross section body considered 
in these comparisons is forebody 2 from Reference 18, and 
the predicted results are from Sample Case 8 in Figure 14 (f ) . 
The comparison -of measured and predicted pressure distribu- 
tions in Figure 31 are in good agreement over the entire 
body surface, even on the windward side. The pressure 
distribution on the body is such that separation does not 
occur over most of the length of the forebody. Separation 
begins at approximately x/L = 0.8, but in contrast to the 
elliptic cross section bodies with their major axis 
horizontal, the vortices are very weak and located near 
the lee side stagnation point. 

Arbitrary Bodies 

The arbitrary noncircular body option in N02VTX has been 
checked out for a number of different cross section shapes to 
verify the numerical transformation. Some of the shapes 
considered are a square cross section with rounded corners, 
a diamond cross section with rounded corners , and a lobed 
or pear shape cross section. The latter cross section is the 
only shape for which pressure data are available, therefore, 
results on this configuration are described as follows. 

The configuration selected for these comparisons is the 
lobed cross section body, forebody 3 from Reference 18, shown 
in Figure 13. The panel arrangement shown in Figure 13(b) 
is used for the predicted results included herein, and the 
specific calculations were made using Sample Case 10. 
Comparison of measured and predicted pressure distributions 
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at two axial stations are shown in Figure 32. At x/L = .14 
in Figure 32 (a) , the agreement is very good on the windward 
side of the body, only fair in the region of minimum 
pressure, and poor near the lee side stagnation point. This 
particular location near the nose of the body has minimal 
vortex interference effects. At the mid-point of the 
forebody length, x/L = 0.5, shown in Figure 32(b) , the agree- 
ment is only fair over the entire circumference even though 
the trend of the data is predicted well. One unusual result 
in Figure 32(b) is the overpredicted pressure on the wind- 
ward side of the body, a situtation which has not occurred 
on any other configuration. The roughness in the predicted 
pressure distribution on the lee side of the body is caused 
by local vortex interference effects. The vortices are 
distributed along the body surface, and they have not started 
to roll up. A sketch of actual predicted vortex positions is 
shown in Figure 33. 


CONCLUSIONS 

An engineering method and the program code NOZVTX to 
predict the vortex shedding from circular and noncircular 
bodies in supersonic flow at angles of attack and roll are 
described. .The coupling of a supersonic panel method, an 
axisymmetric source and doublet model and a two-dimensional 
crossflow vortex shedding method has proved to be a successful 
means to predict the general flow characteristics about 
slender, circular and noncircular bodies at large angles of 
incidence and moderate Mach numbers. Comparisons of 
measured and predicted pressure distributions and associated 
flow fields on a variety of bodies indicate that the 
principal features of the flow phenomena are modeled. The 
predicted vorticity distributions in the wake of a body of 
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revolution in supersonic flow are in reasonable agreement 
with experiment. This, also results in reasonable agreement 
between measured and predicted velocity components in the 
flow field of the body. The ability to model the correct 
flow field near the body leads to the capability to calculate 
wake induced interference effects on fins and other control 
surfaces as well as surface pressure distributions. The 
vortex shedding analysis described herein can be incorporated 
into an overall computation method for missiles and aircraft. 


RECOMMENDATIONS 


In the course of development of the program code NOZVTX, 
several specific areas for needed improvements to the method 
were discovered. These improvements were beyond the scope of 
the present study, but they are noted here as specific 
recommendations for future work. 

The first, and possibly the most important, recommenda- 
tion involves a thorough testing of program NOZVTX. Limited 
numbers of comparisons with experimental results on a 
variety of bodies are included herein, but many other 
comparisons are needed. The program should be tested over 
a wide range of angles of attack and roll and Mach numbers 
to better define the operational limits of the method. The 
program should also be applied to a wide range of noncircular 
body shapes for which data are available. 

A critical weakness in the prediction method appears to 
be the difficulty in calculating accurate pressure 
distributions on the windward side of certain bodies. For 
example, from the comparisons with experiment described in 
a previous section, the method has problems with conical 
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forebodies at all Mach numbers and most shapes at high Mach 
numbers. These difficulties may be caused by nonlinear 
effects induced by the shock wave when it is relatively near 
to the windward body surface. Further investigation into 
this problem is recommended. 

Comparisons of measured and predicted velocity components 
on the lee side of an axisymmetric body at high angle of 
attack indicate the method has a potential capability to 
handle this problem. However, some additional work on the 
details of the vortex cloud are required. This may involve 
a different procedure to blend the discrete vortices when 
velocities in the midst of the cloud are to be calculated. 

With regard to the cloud itself, qualitative results on a 
3:1 missile body are in excellent agreement with experiment; 
but based on one set of experimental results, predicted 
vortex clouds on an axisymmetric body are typically lower 
than those measured. Additional investigation into the 
mechanics of the vortex tracking procedure may explain this 
disagreement . 

The final set of recommendations involve improvements 
to various parts of the program. The Stratford separation 
criteria have proved to work very well even though they are 
based on incompressible flow; however, similar criteria 
or corrections to the existing criteria should be developed 
for compressible flow. In this same area, it has been 
observed in surface flow visualization studies that the 
separation line on a body exhibits characteristics of laminar, 
turbulent, and transitional flow. At the present time, the 
crossflow plane separation may be either laminar or turbulent 
but not transitional or a combination of laminar and turbu- 
lent. It may be possible to use surface flow visualization 
information to correlate the three separation regions on the 
body and build this feature into the prediction method. 
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Figure 3. - Sketch of crossflow plane separation points. 
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Figure 7.- Subroutine Cross Reference Table for Program NOZVTX. 
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Figure 7.~ Continued. 
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Figure 7.~ Continued. 
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Figure 7.- Continued 
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Figure 7.- Continued. 
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Figure 7.- Concluded. 
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Figure 8.- Common Block Cross Reference Table for Program NOZVTX. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Concluded. 
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Figure 9.- Input forms for Program NOZVTX. 
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Figure 9.- Continued. 
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Figure 9.- Continued. 
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Figure 9.- Continued. 
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Figure 9.- Concluded. 
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Figure 14.- Continued. 


121 


1 0 1 

1 

0 

0 1 

1 

0 

0 

0 

0 0 0 

0 

* ! 1 

1 

3 

2 0 







NIELSEN ENGINEERING *- RESEARCH. INC 




PROGRAM iNOZVTX 


SAMPLE Case A 










3: 1 ELLIPTICAL 

CROSS-SECT I ON BODY 



m£F: NASA 

TM-7A079 By 

GRAVES 

MACH NO. = 2.5. 

ALPhAC = 20 DEGREES. 

°H I = *5 

DEGREES 



12.56588 *.P 


16.8 

28.0 


*.0 





?0.0 *s.o 


666*63. 

0.0 


2.5 





P.6 28.0 


1 . * 

1.05 


0.0 


1.0 

O 

O 


.01 0.0 


O 

O 

0.0 


0.0 





30 

0.0 .*20 


.700 

2.100 


3.500 


*.900 

5.600 

6. 300 

7.700 8. *00 


9. 100 

10.500 


1 1.900 


12.600 

13.300 

1*. 700 

\5. *00 16.800 

17.500 

1 8.200 


19.0*0 


19.700 

2u. 300 

21.000 

?2.*00 23.100 

?* . 50n 

25.900 


26.600 


28.000 



o.o .lie 


.196 

.587 


.650 


1.078 

1.162 

1.281 

1 . *65 1.550 


1.631 

1.781 


1.916 


1.977 

2.03* 

2.136 

2.179 2.2S2 


2.279 

2.299 


2.3 09 


2.30* 

2.288 

2.265 

?.2 02 2.165 


2.096 

2.006 


1.970 


1.923 



.260 .280 


.280 

.23* 


.175 


. 158 

• 1*5 

.135 

.128 .119 


.110 

.102 


.093 


. 065 

.076 

.069 

.055 .0*7 


.03* 

.020 


.00* 


-.017 

-.030 

-.0*1 

-.0*7 -.055 


-.057 

-.055 


-.039 


-.03* 



0 . 0 0 .176 


. 29* 

.881 


1.275 


1.61 7 

1.77* 

1.922 

2.19R 2.326 


2 . * * 7 

2.672 


2.8 7* 


2.965 

3. o5 1 

3.203 

3.269 3.378 


3. *18 

3 . ** 8 


3. *6* 


3. -So 

3.*32 

3.397 

3.303 3.2*8 


3.128 

3.009 


2.955 


2.665 



0.00 .059 


.098 

.29* 


.*25 


.539 

.591 

.6* 1 

'.733 .775 


.816 

. 691 


. 958 


• 9ee 

1.017 

1 . 0 68 

1.090 1.126 


1 . 1*0 

1.1*9 


1 .155 


1.152 

1.1** 

1.132 

1.101 1.063 


1.0*3 

1 .003 


. 985 


. 962 



3:i ELLIPTICAL 
l 0 

-? -1 10 30 

CROSS-SECTION RODy 

- 

rfDYbUY 

GEOMETRY INf-uT 











0.0 .*20 . 7013 

2.10 

3. SO *. 

9n 

5.60 

6 

. 30 7. 

70 6 . *0 


9.10 10.50 11 

.90 

12.60 

13.30 1* 

.70 

15. *0 

16. *0 17 

.50 18.20 


19.0*0 19.70 20. 

30 

21.00 

22. *0 23. 

10 

2* . 60 

25 

.90 26. 

60 2e.oo 


0.0 .176 .29* 

.881 

1.275 1. 

617 

1.77* 

1 

.922 2. 

198 2.326 


2 . ** 7 2.672 2. 

87* 

2.965 

3. r >51 3. 

203 

3.269 

3 

.378 3. 

*18 3.**8 


3. *6* 3. *56 3. 

*32 

3.397 

3.303 3. 

2*8 

3.12« 

3 

.009 2. 

955 2.865 


0.0 .059 .098 

.20* 

.*25 

539 

.591 


.6*1 

733 .775 


.016 .091 .958 

.988 

1.017 1. 

0*8 

1 . 09(1 

1 

.126 1. 

1*0 1.1*9 


1.155 1.152 1. 

1 ** 

1.132 

1.101 1. 

063 

1.0*3 

1 

.003 

985 .962 


3:1 elliptical 

CPOSS-SECT ION BODY 

- 

additional 

PANELING InmuT 



0 0-1 

n 0 


(d) Sample Case 4 
Figure 14.- Continued. 
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NIELSEN 

ENGINEERING RESEARCH. I NC 

. 

PROGRAM NOZVTx 


SAMPLE CASE S 







CIRCULAR 

OGIVE 

cylinder 


REE: aFaTL 

TR-70-127 

BY OBEfiKAMPE 


MACH NO. 

= 2.0* 

AlHhac = 1* 

DEGREES. 

RhI * 0 DEGREES 



7.06858 

3.0 

6.0 

39. 

3.0 




15.0 

0.0 

1 760000. 

0.0 

2.0 




1.0 

39. 

1.5 

1.05 

0.0 

1.0 

0.0 


.01 

0.0 

0.0 

0.0 

0.0 




26 








o.O 

.25 

.so 

.75 

1.00 

1.25 

1.50 

1.75 

P. 00 

2. 25 

2.50 

2.7S 

3.00 

3.25 

3.50 

3. 75 

4.00 

A. 25 

A.50 

a. 75 

5.00 

5.25 

5. SO 

5.75 

6. 00 

36.0 







0..0 

.12981 

.25272 

. 36895 

.47870 

.58216 

. b7948 

77082 

. «5630 

.93606 

1 . 01020 

1.07883 

1 . 1 420? 

1 . 1 9Q90 

1.25250 . i 

.29990 

1.34216 

1.37933 

1 .41 146 

1 • 43858 

1.46072 

1.47792 

1.49019 1 

.49755 

1 .500 

1.500 







.a-3333 

.50528 

.47815 

.45185 

.42630 

. 4 0 1 4 5 

.37722 

35355 

• 330 a 1 

.30773 

.2»SA8 

.26361 

.2a209 

. 22i>»9 

.1999o 

1 7 928 

.15883 

.13857 

.11847 

.09851 

.07867 

. 05893 

.03925 

01961 

n.O 

0.0 







21.0 

30. 

3R. 







n. O 1.2 

o. O 1.50 

0.0 l.flO 

0.0 2.10 

0.0 2.*0 

0.0 2.70 

0.0 3.00 

.30 1.2 

.3 1.50 

.3 1 . oo 

.3 2.10 

.3 2. *0 

.3 2.70 

.3 3.00 

. *5 1.05 

.AS 1.10 

.AS 1.2 

.AS 1.30 

. a 5 l.AO 

.as l.SO 

.A5 1.60 

.AS 1.70 

.AS l.SO 

,A5 1.90 

.AS 2.00 

.A5 2.10 

.45 2.20 

.AS 2.30 

.45 2. AO 

.45 2.70 

.AS 3.00 

.60 1.05 

.60 1.10 

.60 1.2 

.60 1.30 

.60 l.AO 

.6 1.50 

.6 1.60 


(e) Sample Case 5 
Figure 14.- Continued. 
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6 

1.70 

6 

1.80 

6 

1.90 

6 

2 . 0-0 

6 

2.10 

6 

2.20 

6 

2.30 

.6 

2.40 

.6 

2.70 

6 

3.00 

75 

1.2 

75 

1.30 

75 

1.4 

75 

1.50 

75 

1.6 

75 

1.7 

75 

1 . 8 " 

75 

1.9 

75 

2.0 

75 

2.10 

75 

2.20 

75 

2.3 

75 

2.40 

75 

2.70 

75 

3.00 

9 

.75 

9 

. H 

9 

. 9 

9 

1.0 

9 

1 . 1 

90 

1.2 

9 

1.3 

9 

1 .4 

9 

1 .50 

9 

1.6 

9 

1 . 7 

9 

1.9 

9 

1.9 

9 

2.0 

9 

2. 10 

9 

2.2 

9 

2.3 

9 

2.40 

.9 

2.7 0 

9 

3.00 

. 05 

. 6 

.05 

. 7 

.05 

. 8 

. 05 

. 9 

.05 

1.0 

.05 

1.1 

.05 

1.2 

.05 

1.3 

.05 

1.4 

.05 

1.50 

.05 

1.6 

. 05 

1.7 

.05 

1 . 80 

.05 

1.9 

.05 

2.0 

.05 

2.10 

.05 

2.2 

.05 

2.3 

.05 

2.40 

. 05 

2.70 

. r 9 

l.nn 


(e) Sample Case 5, Continued 
Figure 14.- Continued. 
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1.20 

0.0 

1.20 

.10 

1 .20 

.2 

1.20 

.3 

1.20 

.4 

1.20 

.5 

1.20 

.6 

1.20 

.7 

1.20 

.8 

1 . 2n 

.Q 

1.20 

1.0 

1.20 

1.1 

1.20 

1.2 

1.20 

1.3 

1.20 

1.4 

1.2 

1 . 5n 

1.2 

1.5 

1.2 

1.7 

1.2 

1 .* 

1 .2 

1.9 

1.2 

2.o 

1.2 

2.10 

i .2 

2.2 

i .2 

2.3 

1.2 

2.40 

1.2 

2.70 

1.2 

3.00 

1 . 3 R 

0.0 

1 . 35 

. 1 

i . 35 

.2 

i . 35 

. 3 

i . 35 

. 4 

i .35 

.5 

' . 3=i 

.6 

i .35 

.7 

’ . 35 

.8 

’ . 35 

. 9 

1 . 35 

1.0 

1 . 35 

1 . 1 

i . 35 

1,2 

1.35 

1 .3 

1 . 35 

1.4 

1.35 

1 . 5n 

1 . 35 

1 .6 

1 . 35 

1.7 

1.35 

1.80 

1 . 35 

1.9 

1.35 

2.0 

1.35 

2.10 

1 .35 

2.2 

1 .35 

2.3 

1 . 35 

2.40 

1.35 

2 . 70 

1.35 

3.00 

1 .50 

0.0 

1 .50 

. 3 

1.50 

.6 

1 .50 

.9 

1.50 

1.2 

1.50 

1 . S0 

1.50 

1.80 

1.50 

2.10 

1.50 

2.40 

1 .50 

2.70 

1.50 

3.00 

51 

? 8.0 


(e) Sample Case 5, Concluded. 
Figure 14.- Continued. 


125 



r 


0 0 0 0 0 0 11 

A 1 0 0 0 0 0 

0 

0 

0 0 

0 0 

0 

0 

NIELSEN ENGINEERING A RESEARCH. ]NC. 
SAMPLE CASE * 



program 

NOZVT* 




CIRCULAR OGIVE CYLINDER REF: NASA TM a- 355 « 3 Y LANDRUM 

MACH NQ . s l.fc. alRmAC * 0 DEGREES* HH I s 0 DEGREES 


7.07 

20.0 

9.0 

20.0 

3.0 




n . 0 

0.0 

=> 00000 . 

0.0 

1.6 




.50 

IP . 5 

0.0 

1.05 

0.0 

1.0 

1.0 


.01 

0.0 

0.0 

0.0 

0.0 




11 









1 . 

2 . 

3 . 

4 . 

5 . 

e>. 

7 . 

a . 

9 . 

2 r • 






0 . 

.32 

.40 

. $4 

1.04 

1.22 

1.34 

1.42 

1 .*<* 

1.50 

1.50 






.343 

.301 

.261 

.221 

.183 

.144 

. 109 

. 0/2 

.0 36 

0.0 

0.0 






51 

2 9.0 

2 ". 







(f) Sample Case 6 
Figure 14.- Continued. 
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1 0 0 1 0 0 1 1 
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0 

NIELSEN ENGINEERING 9 RESEARCH. INC. 
SAMPLE CASE 7 



PROGRAM NOZVTX 




2 : 1 ELLIPTICAL CROSS SECTION FOREBODy REF: NASA TM-80062 BY TOkNSENO 

mach NO.* 1.7. Ai.PHAC * 20. 3- DEGREES. PHI * 0 DEGREES 


13.17679 

14.0 

7.0 

1 4 • 0 

3.60342 




20.36 

0.0 

900855. 

0.0 

1.7 




.50 

14.0 

.50 

1.05 

0.0 

1.0 

0.0 


.05 

39 

0.0 

0.0 

0.0 

0.0 




0.0 

.25 

.90 

. 75 

1.00 

1.25 

1.50 

1. 75 

2.00 

2.25 

2.50 

2.75 

3.00 

3.25 

3.50 

3. 75 

4.00 

4.25 

4.50 

4.75 

5.00 

5.50 

6.00 

6.50 

7.00 

7.50 

8.00 

8.50 

9.00 

9.50 

10.0 

10.50 

I 1 .0 

11.50 

12.0 

12.50 

13.0 

13.50 

1 4 • 0 


0.00000 

0.06764 

0. 13A05 

0.19925 

0.26323 

0.32599 

0.38753 

0.44785 

0.50696 

0.56484 

0.62151 

0.67696 

0.73119 

0. 7<ja20 

0.83599 

0.88657 

0.93592 

0 . 9S40d 

1 .03098 

1 .07668 

1.12116 

1 .20646 

J .28689 

1 . 36245 

1.43313 

1 . 4 9894 

1 .55987 

1.61593 

1.66711 

1.71342 

1 -754B6 

1 .79142 

l.«2310 

1 .84991 

1.87185 

1 .88891 

1.90109 

1.90841 

1.911 


0.27300 

0.26810 

0.26322 

0.25836 

0.25348 

0.24860 

0.24372 

0.23886 

0.23398 

0.22910 

0.22424 

0.21936 

0.2144.8 

0.20960 

0.20474 

0. 1 9986 

0 . 19498 

0.190)2 

0.18524 

0. 18036 

0. 17426 

0.16573 

0.15599 

0 . 1 46?4 

0.1 3649 

0.12674 

0. 1 1690 

0. 1 0724 

0. 09749 

0.06775 

0.07600 

0.0 662*. 

0 . 05849 

0 . 04875 

0.03900 

0 . 02924 

0.01950 

0. 00991 

0.00000 


0.00000 

0.04509 

0.08937 

0.13283 

0.1 7549 

0.21733 

0. 25835 

0.29657 

0.33797 

0.37656 

0.4)434 

0.45131 

0.48746 

0.52280 

0.55733 

0.59105 

0.62395 

0.65604 

0.68732 

0.71779 

0 • 74744 

0.60431 

0.65793 

0.90830 

0.95542 

0 . 99929 

1 .03991 

1.07729 

1.11141 

1 . 1-228 

1 . 16990 

1 . 1942d 

I .21540 

I .23327 

1 . ?4 790 

1 .2SQ27 

1.2674Q 

1.27227 

1 . 27-.00 


0.00000 

0. 0Q018 

0. 1 7874 

0.26567 

0. 35097 

0.4 3465 

0.51671 

0 • 5 ^ 7 1 4 

0.67594 

0.7531 9 

0.82e68 

0.90261 

0.9749? 

1 .04560 

1 . 1 W66 

1 . 18209 

1 .24790 

1 .31208 

1 .37464 

1 .43557 

1.49488 

1 . 60862 

1 . 71586 

1 . 8 1660 

l .9lnS4 

1 .99859 

2.07983 

2.15457 

2.22282 

2. 2S-.56 

2.33981 

2.38855 

2.43080 

2.46655 

?.4958" 

2.51854 

2.53479 

2 . 5--5- 

2.54800 


?:i elliptical c»oss-sect ion 

r Oftt-ODY 

■DY=DY 

GEUME TRY 

INPUT 


u 

-3 -1 18 

29 







0 . 0 

.50 1.0 

1.5 

2.0 2.5 

3.0 

3.5 4.0 

4.5 


5.0 

5.5 6.0 

6.5 

7.0 7.5 

8.0 

8.5 9.0 

9.5 


10.0 

10.5 11. 

11 .s 

12. 12.5 

13. 

13.5 14. 



0.0 

08937 .17549 

.25835 . 

33797 .41434 .48746 

.55733 .62395 .68732 


."'4744 . 

80431 .85793 

.90830 . 

95542 .999291.039911 

.077291.111411.1-228 



1 . 16RO01 . IRA 281 .215401 .233271 .2*790 I . 259271 .267*01. 272271. 27* 

n.O .W87* .35097 .51671 .67594 .82868 .97421 1 . 1 14661 . 2-790 1 . 37464 

1.494 881.608621. 7 1 586 1 . A 1 660 1 . 9 1 084 i . 998592.0 798 32. 1 5*5 72 . 22282? . 28456 
2.339812.388552.430802.466552.495802.518542.534792.544542.548 

? : 1 ELLIPTICAL CROSS-SECTION FOREtiODY - ADDITIONAL PANELING INPUT 
9 0-1 

o 11 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 8.0 10.0 12. 

I*. 


(g) Sample Case 7 
Figure 14.- Continued. 
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1 0 0 1 0 0 1 1 

* 1 1 1 3 2 0 

0 

0 

0 0 

0 0 

0 

0 

NIELSEN ENGINEERING 6 RESEARCH. INC. 
Sample CASE 8 



PROGRAM 

NOZVTX 




2:1 ELLIPTICAL CROSS SECTION FOREBOOY REF : NASA TM-80062 6* TOWNSEND 


MACH NO 

.= 1 . 7 . 

ALPhac * 20 . 

36 DEGREES 

* PHI = 0 

DEGREES 



13.17679 

1«.0 

7.0 

1*.0 

3 . 603*2 




20.36 

0.0 

900855 . 

0.0 

1.7 




.50 

1 *. 0 

.50 

1-05 

0.0 

1.0 

0.0 


.05 

o ts 

0.0 

0.0 

0.0 

0.0 




37 

n.o 

.25 

.80 

.75 

1.00 

1.25 

1.50 

1.75 

p . 00 

2.25 

2.50 

2.75 

3.00 

3.25 

3.50 

3.75 

4.00 

*.25 

*.80 

*.75 

5.00 

5.50 

6.00 

6.60 

7.00 

7.50 

8.00 

8.50 

9.00 

9.50 

10. 0 

10.50 

11.0 

11.50 

. 12.0 

12.50 

13.0 

13.50 

1*.0 


0.00000 

0 . 0676 * 

0 . 13*05 

0.19925 

0.26323 

0.32599 

0.38753 

0.**785 

0.50696 

0 • 56 * 8 * 

0.62151 

0.67696 

0.73119 

0 . 78*20 

0.83599 

0.80657 

0.93592 

0 . 98*06 

1.03098 

1 . 07668 

1.12116 

1 . 206-6 

1 . 28689 

] . 362*5 

1.43313 

1 .* 989 * 

l . 559 B 7 

1.61593 

1.66711 

1 . 713*2 

1 . 75*06 

1 . 791*2 

1.82310 

1 . 8*991 

1.87185 

1 . 88891 

1.90109 

1 . 9 G «* 1 

1.911 


0.27300 

0.26810 

0.26322 

0.25836 

0 . 253*8 

0 . 2*860 

0 . 2*372 

0.23666 

0.23398 

0.22910 

0 . 22 * 2 * 

0.21936 

0 . 21**8 

0.20960 

0 . 20 * 7 * 

0. 19906 

0 . 19*98 

0.19012 

0 . ] 852 * 

n . 16036 

0 . 1 7*26 

0 . 16573 

0. 15599 

0 . 1 * 62 * 

0 . 1 36*9 

0 . 1267 * 

O . 1 1699 

0 . 1072 * 

0 . 097*9 

O .08775 

0.0 7800 

0 • 0 b « 2 * 

0 . 058*9 

0.0*875 

0 .03900 

0 . 0292 * 

0.01950 

0.00991 

0.00000 


0.00000 

0 . 0*509 

0.08937 

o . 13283 

0 . 175*9 

0.21733 

0.25035 

0.29857 

0.33797 

0.37656 

0 . 4 1 * 3 * 

0.*5131 

O .* 07*6 

0.52280 

0.55733 

0.59105 

0.62395 

0 . 6560 * 

0 . 68 73 p 

0.71779 

0 . 7 * 7 ** 

0 . 60*31 

0.85793 

0.90830 

0 . 955*2 

0.99929 

1.83991 

l .07729 

1 . 111*1 

1 . 1*228 

1 .16990 

1 . 19*20 

1 . 215*0 

1.23327 

1 . 2 * 790 

1.25927 

1 . 267*0 

1 .27227 

1 . 27*00 


0.00000 

0.09018 

0 . 1 787 * 

0.26567 

0.35097 

0 . * 3*65 

0.51671 

0 . 5971 * 

0 . 6759 * 

0.75319 

0.92868 

0.90261 

0 . 97*92 

1 . 0*560 

1 . 11*66 

1 .10209 

1 . 2*790 

1 .31208 

1 . 37 * 6 * 

1.43557 

1 .* 9*88 

l . 60*62 

1.71506 

1 .81660 

1 . 9108 * 

1,99859 

2.07983 

2 . 15*57 

2.22282 

2 . 28*56 

2.33901 

2 . 3 B 655 

2. *3080 

2.46655 

2. *9580 

2 . 5185 * 

2 . 53*79 

2 . 5 * * 5 * 

2 . 5*600 



? : 1 ELLIPTICAL CROSS-SECTION fq*£p 0D Y - wDY-OY GEOMETRY INPuT 
0 0 

-3 -1 18 29 ■ 


0.0 

.50 

1 . o 

1 .5 

2.0 

2.5 3.0 3.5 

*.0 *.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 6.0 6.5 

9.0 9.5 

10.0 

10.5 

1 1 . 

11.5 

12 . 

12.5 13 . 13.5 

1 *. 

0.0 

.08937 

. 175*9 

.25835 

.33797 

.* 1 * 3 * .* 87*6 .55733 

.62395 .66732 

. 7 * 7 ** 

. 00*31 

.05793 

.90830 

. 955*2 

. 999291 . 039911 .077291 

. 111 * 11 . 1*228 


1 . 169901 . 1 9 * 28 1 . 2 1 540 I . 23327 I . 2 fc 790 1 . 259271 . 267 * 0 1 . 2 7227 1.274 

0.0 . 1787 * .35097 .51671 . 6759 * .82868 . 97 * 2 1 1 . 1 1 . 66 1 . 2 * 790 1 . 3 7 * 6 * 

1 .* 9 * « 81 . 608621 . 715861 . 816601 . 91 09 * 1 . 998592 . 079832 . 15 * 572 . 222822 . 28*56 
2 . 3399 12 . 398552 . * 3 08 02 . * 66552 . * 958 02 . 5 195 * 2 . 53 * 782 . 54 * 5 * 2 . 5* 8 

?:i elliptical cross-sect»on forehody - additional panel I ng input 
0 0-1 
- 1 « 11 

0.0 5 . 10 . 18 . 26 . 30.5 * 5 . 58.75 72.5 107.5 

121.25 135 . 1**.5 15 *. 162 . 170 . 175 . 160 . 

0.0 1.0 2.0 3.0 *.0 5.0 6.0 6.0 10.0 12 . 

1 *. 


(h) Sample Case 8 
Figure 14.- Continued. 
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1 0 0 

1 0 2 

1 1 

0 

0 0 0 

0 0 

4 1 1 

13 2 

0 




SAMPLE CASE 9 
3: 1 ELLIPTICAL 

CROSS SECTION 

CONE 

«£F: 

NASA T M- 78808 

8Y TOWNSEND 

mach no.* 1.7. 

ALPMAC * 20.35 

DEGREES. 

PHI * 

0 DEGREES 


27 . 1896 6.79* 

7.0 

l*.o 

6.794 



?0 . 35 0.0 

1130000. 

.25 

1.70 



.50 7.0 

0.0 1 

.05 

0.0 

1.0 

.50 

.01 0.0 
2 

0.0 1*.0 

0.0 3.397 

.2*27 .2*27 

n.O 5.0055 

0.0 1.69R5 

0.0 

0.0 

0.0 




3 : 1 CONICAL ELLIPSE - WDY8DY GEOMETRY INPUT 
n 0 

-3 -1 IP 2 

0 . 1 4. 0 

o » Q 5.0 955 
0.0 1.6985 

3:1 CONICAL ELLIPSE - ADDITIONAL PANELING INPUT 
0 n~l 
-1 A 5 

0.0 13,85 27. 69 *1.54 55.39 69. S» 77.00 81.50 97.00 93.00 

9a . 50 103. 110.5 124.62 13P.46 152.31 166.15 190. 

0.0 .50 6.25 12.0 1 4 . 0 


(i) Sample Case 9 
Figure 14.- Continued. 
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2 

0 

0 l 

0 2 

I 1 

0 0 

0 

0 0 '0 

0 

4 

1 

1 1 

3 2 

0 





NIELSEN 

ENGINEERING k RESEARCH* INC. 


PROGRAM NOZVTX 


SAMPLE 

CASE 10 







L08ED 

CROSS SECTION FORE6QDY 

REP: NASA 

TM-80062 

bY TOWNSEND 

m*ch no 

.* 1.7. ALPHAC = 20. 

36 DEGREES 

. PHI * 0 

DEGREES 



15.05941 


14.0 

7.0 

1 4 • 0 

4.37884 




20.36 


0.0 

72876*. 

.075 

1.7 




.50 


7.0 

0.0 

1.05 

0.0 

2.0 

1.0 


.05 


0.0 

0.0 

0.0 

0.0 




29 









0.0 


.50 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 


4.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 

8.0 


8.5 

9.0 

9.5 

10. 

10.5 

11.0 

11.5 

12.0 


12.5 

13.0 

13.5 

1*.0 




0.0 


.15359 . 

.38160 

.44403 

.58087 

.71212 

.63779 

.95787 

1.07237 


1.18128 

1 . 28461 

1.36235 

1.47451 

1.56108 

1.64207 

1.71747 

1 . 78728 


1.85151 

1.91016 

1.96322 

2.01069 

2. 05258 

2.08889 

2.11961 

2. 14474 


2.16429 

2. 1 7826 

2.16663 

2. 1 8942 




.31 


.30160 

. 29044 

.27927 

.26809 

.25692 

.24575 

.23458 

.223*1 


.2122* 

.20107 

. 18990 

.17873 

.16756 

. 15639 

. 1*521 

.13*04 


.12288 

.11171 

. 1 0053 

.08936 

.07820 

.06703 

. 05585 

. 04468 


.03351 

.02234 

.011J7 

0.0 





20 1 

50 

2 * 

.00000 -.17874 

.01558 -.17809 

.03105 -.1761? 

.04629 -.17275 

.06110 -.16786 

.07524 -.16134 

.08840 -.15312 

.10026 -.14319 

.11049 -.13168 

.11882 -.11882 
.12508 -.10*96 

.12926 -.09051 

.131*7 -.07590 

.13194 -.06152 

.13097 -.04767 

.12891 -.r.3454 


. 12607 
. 11*00 
.10000 
.08365 
.07647 
.07212 
.06431 
.05476 
.04333 
.03011 
.01546 
•0000O 


-.02223 
.01900 
.05800 
.09970 
. 11206 
.12492 
.13791 
. 15044 
.16171 
.17077 
. 17666 
.17874 


LOBED 
n 1 
1 1 

CROSS SECTION 
20 29 

POOY - 

wOYBQY GEOMETRY 

INPUT 



n. 

.5 

1.0 

1.5 

2. 0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

S .5 

6.0 

6.5 

7.0 

7.5 

8.0 

8.5 

9.0 

9.5 

10. 
0.0 
n . n 

10.5 11.0 

0.0 
n.n 

11.5 

12. 

12.5 

13. 

13.5 

14. 



0 


(j) Sample Case 10 
Figure 14.- Continued 



o.O 
0.0 
0.0 
0.0 
0.0 
0.0 
o.O 
0.0 
0.0 
0.0 
o.O 
0.0 
0.0 
0.0 
0.0 
o.O 
0 . 0 
o.O 

.00000 
.01550 
.03105 
. 0 * b 29 
.05110 
. 056*0 
. 110*5 
•12508 
. 131*7 
. 13097 
.12607 
. 11*00 
. 1 0000 
.08365 
• 07212 
. 05*76 
. 04333 
.03011 
. 015*6 
.00000 
.00000 
. 03059 
.06098 
.09069 
.11957 
. 17359 
.21696 
. 2*561 
.25816 
.25719 
. 2*756 
.2239 
. 196 * 

. 16*27 
. 1*162 
. 10752 
.08509 
.05913 
.03035 
0.00000 
. 00000 
. 0*504 
.08970 
• 13381 
.17662 
.25557 
.31942 
.16160 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

17674 
-.1 7809 
-.17612 
1 7275 
-. 16786 
-. l c 312 
-.13168 
-. 10 ** 96 
-.07590 
-.04767 
-.02223 
.019 
, 05 « 

. 09970 
. 12*92 
. 150 ** 
.16171 
.17077 
. 17668 
. 1 787 - 
-.35097 
. -. 3*970 
- . 14584 
-.33922 
-.32962 
-.30067 
-.25857 
-.20609 
-. 1*905 
-.09361 
- . 04 36 S 
.0373 
.1139 
. 19577 
. 2*530 
. 295*1 
.11755 
. 3353 * 
. 3*693 
.35097 
-.51671 
- . 51 * 8 * 
-.50915 
-.49940 
-.48527 
-. 4*265 
-.38067 
-. 303*2 


(j) Sample Case 10, Continued. 
Figure 14.- Continued. 
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.38006 
.37863 
. 36*^7 
. 12956 
. 28909 

. 2 * 18 * 

.20050 

.15830 

.12527 

.08705 

. 0**69 
0.00000 
.00000 
.05892 
. 117 ** 
. 1 7505 
.23106 
. 33*32 
.41786 
. 4730 * 
.49710 
.49532 
.47678 
.-3113 
. 37818 
.31637 
.27276 

.20708 
.16387 
. i 138 *> 
. 05846 
0.00000 
.00000 
.07224 
. 14398 
. 21*61 
.23327 

.40987 

.*1223 
.57992 
.60953 
.60724 
.50452 
.52856 
. 4636 * 
.38785 
. 33*39 
.25387 
.20090 
•13961 
.07167 
0.00000 
. 00000 
. 09*99 
.16939 
. 252*8 

.33326 

.*8220 

.60268 

.68226 

.71710 

. 714*0 

.68767 

.*2182 

.54545 

.45630 


-. 219*3 
-.13781 
-. 06*27 
. 05*93 
. 167*7 

.28821 

.36114 

. 43*92 

.*6750 

.49369 

• S 1076 

.51671 

-.67594 

-.67350 

-.66606 

-.65331 

-.634f»2 

-.57907 
-.49798 
-.39692 
-.28705 
-.18028 
- . 08407 
.07185 
.21935 

. 37703 
.47243 
.56895 
.61158 
.64503 
.66817 
.67594 
-.82668 
-.82569 
-.81656 
-.00093 
-.77827 
-.70991 
-.61051 
- .48661 
-.35101 
-.22 102 
-. 1 0307 
.08809 
.26891 

. *6222 
. 5791 a 
.69751 
. 74977 
.79177 
.81915 
.82868 
-.97492 
-.97139 
-.06066 
-. 9*227 
-.91561 
-.83519 
-.71825 
-. 572*9 
-.41402 
-.26002 
-.12125 
. 10364 
31636 
.54 373 


(j) Sample Case 10, 


Continued 


Figure 14.- Continued 


•39340 

.66139 

*2966 7 

.62060 

*2363*5 

.86206 

. 16425 

.93149 

. 68*31 

.96371 

0.00000 

.97492 

.00000 

-1.11466 

.00717 

-1.11063 

.19367 

-1.09835 

.26667 

-1.07733 

.36102 

-1.0*685 

.55131 

-.95491 

.68906 

-.82120 

.78005 

-.65*54 

.81968 

-.47338 


-.29729 

.78624 

-.13863 

. 7 1005 

. 1 1849 

•62364 

.36171 

.5217ft 

.62174 

. 44979 

.77906 

•34J48 

.93822 

.27023 

1.00851 

•18779 

1.06500 

. 09640 

1.10104 

0.00000 

1.1U66 

.00000 

-1.2*790 

. 1 0 87o 

-1.24339 

.21682 

- 1.22964 

.32317 

-1.20610 

.42657 

-1.17198 

.61722 

-I . 0690S 

. 77143 

-.91936 

. 8733ft 

-.73278 

.91 769 

-.52994 

.91 444 

-.332*3 

.8802? 

-.15521 

• 7 9S93 

. 1 3265 

.49818 

. *0*95 

• 584 08 

.69608 

.50355 

.07218 

. 39230 

1.05036 

. 30253 

1. 12906 

.2102* 

1. 19231 

•10792 

1.23355 

0.00000 

1 .24790 

. 00000 

- 1 . 37464 

•11983 

-1 . 36967 

.23684 

-1.35453 

. 35600 

-1.32060 

. 46969 

-1.29101 

.67990 

-1.17762 

.64976 

-1.01273 

.96199 

-.60721 

l.onn 

-.58376 

1.00731 

-.36663 

.96961 

-.1 7097 

' q 76 76 

14613 

76909 

4*607 

.64336 

. 76675 

.55*70 

.96076 

.42113 

1. 1S70* 

.33326 

1.2*373 

.23159 

1 .31340 

. 11068 

1 • 35883 

n.ftftoftft 

J . 37464 


(j) Sample Case 10, Continued. 
Figure 14.- Continued. 
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.00000 
.13031 
.25973 
. 3871 * 
.51090 
.73937 
. 92*11 
1 . 0 * 61 * 
1 .09955 
1 . 095*2 
1 . 05**3 
. 953*5 
'.*3636 

.69965 
.60322 
.*5796 
. 362*1 
. 2518 * 
.12925 
0.00000 
.00000 
. 1*023 
. 279*9 
.*1659 
. 5*987 
.79563 
. 99**2 
1.12573 
1.18321 
1.17877 
1 . 13*66 
1.02600 
.90000 

.75280 

. 6*911 

.*0281 

.38998 

.27101 

.13912 

0.00000 

.00000 

. 1*958 

.29813 

.***37 

.58653 

. 8*867 

1 . 0607 ? 

1 .20078 
1.26210 
1.25735 
1.21030 
1.09*00 
.96000 
,e0308 
.69239 
.52566 
.*1598 
.28907 
. i * e 39 
0.00000 
.00000 
.15836 
.31563 
.*70*5 
.62096 
. soften 


- I . * 9*88 
-l.*89*7 
- 1.47301 
- 1 • -***8 1 
— i . *0 393 
- 1.28063 
- 1.10131 
-.87781 
-. 63*83 
-.39870 
-. 1859 ? 
.15891 
.*8509 
. 8338 ? 
1 . 0**80 
1 .25825 
1.35252 
1 .42828 
1 .*7768 
1 .*9488 
- 1.60862 
- 1.60280 
- 1.58509 
- 1 . 55 * 7 * 
- 1.51076 
- 1 . 37 B 07 
- 1.18511 
- . 9**60 
-.6831 3 
- . * 290 * 
-.20007 
.17100 
.52200 
.89726 
l . 12*30 
1.35398 
1 .* 55*3 
1.53696 
1.59012 
1.60862 
- 1.71586 

- 1.70965 

- 1.69076 
- 1.65839 
- 1 . 611*7 
- 1 .* 699 * 
- 1 . 26*11 
- 1.00758 
-.72867 
“ . * 576 * 
-. 213*1 
. 182*0 
.55680 
.95708 
1 .19925 
1 .***25 
1 . 552*6 
1 . 639*2 
1.69612 
1.71586 
- 1 .81660 
- 1.81003 
- 1.79003 
- 1.75576 
- 1 . 70609 
- 1 . 556?4 


(j) Sample Case 10, Continued 
Figure 14.- Continued. 



1.12299 

1 .2712* 
1.33620 
1.33117 
1.26136 
1.15665 
1.01636 
• 65023 
.73304 
.55653 
.44040 
. 30605 
.15710 
0.00000 
.00000 
.16657 
. 33200 
.49486 
.65318 
.94511 
1.18125 
1.33723 
1.40551 
1.40023 
1.34783 
1.21676 
1.06909 
. 69434 
.77107 
.58540 
.46325 
.32192 
.16525 
0.00000 
.00000 
.17422 
. 34725 
-51 758 
.68317 
. 9885] 

1.23540 
1 . 39864 
1.47005 
1 . 46*53 
1.40972 
1.27473 
1.11618 
.93541 
.80647 
.61228 
.48452 
. 33671 
. 17284 
0.00000 
.00000 
.18130 
.36137 
.53862 
.71094 
1.02869 
1.28572 
1.45540 
1.52981 
1.52406 
1.46703 

1 . -»?<sss 


-1.33833 
-1.06673 
-.77U5 
-.48451 
-.22594 
.19311 
• 5B949 
1.01327 
1.26966 
1.52904 
1.64361 
1.73567 
1.79571 
1.81660 
-1.91064 
-1.90393 
-1.88289 
-1.84685 
-1 . 79459 
-1.63698 

- 1 .40776 
- 1.12207 

-.81147 
-.50964 
-.23766 
.20313 
.62007 
1 . 06583 
1 . 33553 
1.60837 
1 . 72888 
1.82572 
l . 88887 
1.91084 

- 1 . 99889 
-1.99136 

- 1 . 96935 
-1 .93165 
-1.87700 
-1.71216 
-1.47241 
-1.17360 

- . 84674 
-.53304 
-.24857 
.21245 
.64855 
1.11478 
1 . 39685 
1.68222 
1.80826 
1.90955 
1.97560 
1.99659 
-2.07963 
-2.07231 
-2. 04940 
-2.01017 
-1 .9S330 
-1.78175 
-1.53226 
-1.22130 
-.86324 
-.55471 
-.25868 

. ??l 09 


(j) Sample Case 10, Continued 
Figure 14.- Continued. 


1 . 16364 
.97343 

.83926 
.63717 
.50422 
.35039 
.17987 
0.00000 
.00000 
.18782 
.37435 
.55798 
.73649 
1 . 0656ft 
1.33192 
1.50780 
1 . 50*79 
1 .57883 
1.51975 
1 . 37*22 
1 .20545 
1 . 008*1 
. 869*2 
.66007 
. 5223 * 
.36299 
. 1 8633 
0.00000 
.00000 
.19377 
.38621 
.57565 
.75982 
1 . 099*1 
1 . 37*11 
1.55556 
1 . 63*99 
1 . 6288 * 
1.56789 
1.41775 
1.24364 
1 . 0*036 
. 89696 
.60097 
.53888 
. 37**8 
•19223 
0.00000 
.00000 
.19915 
. 3969 * 
.59165 
. 78093 
1 .12995 
1 .*1228 
1.59877 
1 . 680*0 
1.67409 
1 .61144 
1.45713 
1.27818 
1.06925 

.92167 
.69989 
.55365 
. m* a n 


.67491 
1.16009 
1.45364 
1.75061 
1.88177 
1.96718 
2.05591 
2.07983 
- 2 . 15*57 
- 2 . 1467B 
- 2.12305 
- 2 . 082*1 
- 2.02350 
- 1.84578 
- 1.58732 
- 1.26519 
-. 91*98 
-. 57*65 
-.26797 
. 2290 * 
.69916 
1.20178 
1 . 5O5«0 
1.81352 
1 . 9 * 9*n 
2.05859 
2. 1 2979 
2 . 1 5 * 5 ? 
- 2.22202 
- 2 . 21*78 
- 2.19030 
- 2 . 1*837 
- 2 . 0B759 
- 1 . 90 * 2 * 
- 1 .63760 
- 1.30527 
-. 9*396 
-.59285 
-. 276*6 
.23629 
.72131 
1.23980 
3 .55357 
1 . 9709 * 
2.0U1* 
2.12379 
2.19725 
2.22202 
- 2 . 28*56 
- 2.27630 
- 2 . 2511 * 
- 2.20805 
- 2 . 1*558 
- 1 . 9571 * 
- 1.68309 
- 1 . 3 * 1S3 
-.97018 
-.60932 
-. 28 * 1 * 
. 2*285 
. 7*135 
1.27429 
1.59673 
1.92293 
2.06701 

? . 10270 


(j) Sample Case 10, 


Continued 


Figure 14.- 


Continued. 


.19757 
0.00000 
. 00000 
.20397 
. *0654 
.60595 
.79981 
1.15728 

1 .44643 
1 .63743 
1.72104 
1 . 71457 
1.65041 
1.49236 
1.30909 
1.0951 1 
.04416 
.71601 
.56725 
.39419 
. 20236 
0.00000 
. 00000 
.20622 
.41501 
.61650 
.61647 
1.18139 
1 .47657 
1.67154 
1 . 75600 
1.75029 
1 .60479 
1 .52345 
1 . 3363 b 
1.11793 
.96303 
.73175 
.67906 
. 40240 
.20657 
0.00000 
. 00000 
.21190 
.42235 
.62952 
.83091 
1 .20228 
1.60268 
1.70111 
1.70797 
1 . 70120 
1.71459 
1.55040 
1 . 36000 
1 . 1 3770 
.90000 
.74469 
.58931 
.40952 
.21022 
0.00000 
.00000 
.21501 
.42056 
.63077 


2.25029 
2.20456 
- 2.33901 
- 2.33135 
- 2.30550 
- 2.26145 
- 2.19746 
- 2.00447 
- 1.72379 
- 1.37397 
-.99364 
-.62405 
-.29101 
.24073 
.75927 
1.30510 
1.63534 
1.96943 
2.11699 
2.23557 
2 . 312^0 
2 . 339 fl 1 
- 2.38656 
- 2. 37992 
- 2.35361 
- 2.30856 
- 2.24324 
- 2.04623 
- 1 . 75970 
- 1.40259 
- 1.01434 
-.63705 
-.29707 
. 253«1 
. 77509 
1.33229 
1.66941 
2.01046 
2.16110 
2.28215 
2.36108 
2.38055 
- 2.43080 
- 2.42201 
- 2.39524 
- 2 . 34 S 39 
- 2.20292 
- 2. 08242 
- 1.79083 
- 1.42740 
- 1.03228 
-.64032 
-.30233 
.25840 
. 78880 
1.35586 
1.69894 
2.04602 
2.19932 
2.32251 
2.40284 
2.43080 
- 2.46655 
- 2.45763 
- 2.43046 
- 2. 30394 


(j) Sample Case 10, Continued 
Figure 14.- Continued. 



. 8*313 
1.21996 
1 . 52*76 
1.72612 
1 . 61*26 
1 . 607 ** 
1.73961 
1.57320 
1.38000 
1 . 15**3 
.99531 
• 7556 * 
.59797 
.* 155 * 
.21331 
0.00000 
.00000 
.21756 
.* 336 * 
. 6*635 
.65313 
1 . 23**3 
1 . 5*286 
1 . 7*659 
1.83577 
1.82887 
1 . 760 ** 
1.59195 
1.39636 
1.16812 
1.00711 
. 76*60 
.60506 
.* 20*7 
. 2158 * 
0.00000 
.00000 
.21955 
.*3759 
. 6522 * 
.86091 
1 . 2*568 
1.55692 
1.76251 
1.95251 
1 . 8 * 55 * 

1 . 776*8 
1,60636 
1 .*0909 
1.17877 
1.01629 
.77157 
.61058 
.* 2*30 
.21781 
0.00000 
.00000 
.22096 
.** 0*1 
. 656*5 
. 666*6 
1.25372 
1.56697 
1.77388 
1 . 86**6 
1 .887** 


- 2 . 316*9 
- 2 . 1130 * 
- 1.81716 
- 1 .**839 
- 1 . 0 * 7*6 
-.65786 
-.30677 
.26220 
. 800*0 
1.37588 
1.72392 
2.07611 
2.23166 
2.35667 
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(j) Sample Case 10, Continued. 
Figure 14.- Continued. 
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(j) Sample Case 10, Concluded. 
Figure 14.- Continued. 
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Figure 15.- Output for sample case 
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Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15. - Continued. 
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Figure 15. - Continued. 
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Figure 15.- Continued. 



BODY P«NFL Cf.MTnOIO POlMl rf)l>P|> 1 1 'l A T f 


00 

Q) 

tr> 

fC 

* 


JC 


> J* * x. z: C ^ J' T O' St X) o c 

— <* c* — c. — r c «— cr — a • b ri a- m © m rv rv cc cr *n <V *n e <n fv »*■. ® tT^c — iTCiT — >c lT it c « if ■ c if ■ ^r — 

c* x _o >jr o c x cr ?i^iooirc5'0'Oflf f icwncrO'iMnnof f incj'-N4-CH4K--N.vj — « c ^ ^ 
a. r-'r^^ ^c: r ^ H r r\j o (\i r r. r n n n it* © lP ^ ^ f^r^r-Lrcirf^f^p'-— a c r- ^ c C <j <r O' z: *— >c 

o ©©ceccccr©©©©© occoooooooooocooocooooo'Mno— ^ — — i\i <\j f\j rvj — ^ © «-■ «— ?\j ‘V 


v, i t i i i i i i t i i i i i » i « i t r i i i l i 


— o x ^ x ■* c> — — « a* 

— — n-^ipcmc m o a lT c ^n (? — nc 

*■ iTr\^* — '***-ir -#'X:a^rv^r\*£ <* x >**crx©**xx* < c f\ r * * <r ;r 

a. o'\.^xex-<rr\.©er\j^x©cc — >» c >i •* ix ^ 

*_* © c © © ~ c c cc©cec~e©e©e©©~-r\j«-«c©ce©©~*cvj— ©©©o-*rvj.#i/',r(\,^-o©.— rv^iT<a<\ — ec 


> i i i i i i i i t i i i i l i i i l i « i i i i i i t 


©ccoeecco©eecc©©e©ff'^--*tne , if'*’^^t>*— 

cccceccoccoocoooooift«o^N^N^^if. i^o^soN-^'CinccCiniTto^^^^'r^inirir^ c ^lt 

COCCO©©©©©©©©©©©©©^-* — — ^ ^ 

x xxx xxx x axxxxxxxff©x^^f^^» s *^^»^r^^ , **^^^^^-^f^^'"' — ~ ^ * 

u rv rv, r\, rv <v rv rv r\. \ \ \iv fv <v r\j rv iv a. iriinirininiriniPiPiniriririnirinirinkririinirinir LTiT'iriiPif'iPLpiPir iTlt irx 



2 * — \r' 4 ircNxO'C«(\.-fn<irc^x^©^»\jn^ircNCcro^(Vjr) 4 ir^Kxcrc — Mr'’#LT^)Nfic^o^^r 4 if 

•_ ^#.^*^^^p,^. < -^f^fv»fVjr\/(\.f\ifv ; rofv,fyf r )r‘mrmnrr)(nn^ 4444 ^ # 4 ^ '•iTiririririr 

c 

a 


148 


Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15. ~ Continued. 
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Figure 15.- Continued. 


THf ITERATION CONVtROFD AFTfR 7 J TF R A I 1 OnS WITH A TtSf CRUFRIUN RF .0010000 


UJ 

> 

o 

«n 


© o 

© no 

<c 

£ 

m 

o 

4 



© 



























in 


r- 


in nj 

.n 

o 

f\J 

ru 

rn 

B 



























O 


m 

O' 

£ 

o 

in 

m 

r- 

o 


n- 



























©> 

O 


h- 

© ru 

m 

o 

nj 


h- 

« 



























©i r- 

m 

-C 

o 

© 

© 

® 

o 

c> 

O 

© 



























1 


• 



1 


i 

— 

1 

— 

1 



























© 

O' 

<r 

© 


© 

© 

n 

m 

f*- 


« 



























© 

— t 

© O 

rv r«- 

rv 


»-* 

rg O' 

o 



























o 

4 

♦ 

© 

© 

rvj 


o 

m 

o 

O 

© 



























in 

O 

4 


o 

■4 




m 

r- 

in 



























© 

© 

n 


n 


r*- 


o 

© 

c 

o 



























f 


1 


i 



l 


i 


1 




-J 






































4 


no 

■n 

'U 


£> 

X 



r* 

n 

4 


X 

© 

p— 

'V 

nr, 

ru 

£5 

4 

© 
















X 


m 

4 


>C 

fs. 

m 

© 

ru 

O 

o 

ru 

© 



© 


4 




V 


















o 

o 

4 


X 


iT 

■n 



iT. 

n' 


n 



© 

O 

X 



m 


© 

•— « 

© 

© 

© 

r 

X 

X 

r- 

© 




c 


sD 



X 

ru 

4 

O 

© 

e 

© 

© 

o 

4 

ru 

X 

f^. 




© 


m 

4 


© 

O 

ru 


in 


4 


in 




z 


4 

4 

4 





ru 

'u 

ru 

nj 




4 



4 

4 

4 

4 


*U 

4 

r- 

“n 

r*. 

“VJ 

-n 

4 

'n 

o 

« 



























4 


©> 

O 


p*-. 

4- 

rv 

4 

no 

•n in 











1 

1 

i 

i 

i 


1 

1 









© 

1 

© 

in 

i 


© 

i 

p*>* 

“ 

r- 

i 

o 

X 

i 

o 

o 

a 



























© 

in 

lT 

ru 


rn 

£ 


o 

© 

ru 

© 



























4 

©> 


X 

vi © 

'VI 

x rvj 

X 

© 











© 







© 









© 

l n 

4 

4 

4 

in 

c 

rv 

© 

-0 


© 




_J 

> 






» 







1 









mm 

o 

© 

o 

■4- 

ru 


o 

© 

>D 





< 

►— 






UJ 







UJ 









© 


© 

© 

© 

X 



X 

X 

O 

o 





»— * 

4 

O 


p*- 

© 

m 

4 

r^ 



r^ 

4 

") 

© 



o 

4 



o 


• 




• 

* 

• 

• 








VJ 

O 

4 

r- 

© 

© 

in 

o 

iT 

o 

o 

in 

o 

s ■ 


e 

r^. 

4 

o 


•w 

— 

i 


i 


i 


1 



1 


1 




♦— 

£ 

f\J 

O 

X 

o 

© 

Tv 

o 

iT 

c: 

c* 

s 

©■ 

rv 

© 



y 

ru 

V 

x 

© 
















a: 

-1 

in 

4 

*n 

X 

o 

O 


O 

© 

c 

o 


O 

y 

© 

*n 

4 

© 

n 

4 

'n 

















u. 1 

4 

4 

4 



o 



ru- 

V 





rn 


4 

4 


4 

















> 

> 

9 





• 










• 



■ 



© 


X 

o 

—• 

i"*, 

o 

o 

r»- 

— * 

J' 







1 

1 

1 

1 


I 

1 

1 

1 

t 

k 

i 

1 


1 

l 

1 

« 

i 

1 

i 

© 


X 


a 

© 

•4 

*v 

© 

iT 

O 

o 

o 


























4 

ru 

4 

4 

© 

X 


rv 

h- 

© 


o 


























o 


4 

o 

© iT 

n> 

o 

*c 

O 

an uT- 

p*> 


























T 

**“ 

© 

X 

in c 

c 

© 



o 

o 

c 


























t 




i 


i 



1 


i 

"• 












_ 









_ 






















o 

© 





© 

© 

© 

- 

© 

© 





© 


~ 

© 


ru 

o 

i 

O 

— 

— 


X 

X 

o 

© 

© 

m 



— i 


UJ 

UJ 







Lku 

Jj 







uJ 

uJ 




4 

© 

r- 

m 

o 

4 

rn 

r- 

n 

r- 

© 

X 

n 



4 


'■n 

X 

4 

f- 


4 

ru 




o 

ru 

4 


r- 

4 

X 

-n 


4 

4 

© 

J* 

rn 

P- 

rv 

o 

■4 

rr 

© 

m 





a 

u 

4 

lT 


4 

r - , 

r- 

r- 

© 

4 

4 

© 



r-. 

4 


© 

4 

c 

© 


© 



o 

^n 

X 

© 

4 

no 

c 

iT 

© 




UJ 

o 

o 


P^. 

ru 

•■u 

4 


4 


f- 

4 

rn 

4 

ru 

ru 



0 

X 

V- 


o 

©> 

© 


© 

X 

x> 

O 



O 

© 

© 


© 


-J 

4 



4 

© 

rn 



4 

4 

f*-< 


-n 

c 

4 

© 



4 

















© 

4 

uJ 

-U 

X 


m 

rn 


o 

4 



4 

o 



-n 


X 

'U 

V 

r 

— » 



1 


i 


1 


t 



l 



© 

-J 

> 

# 

• 

# 

• 

• 

1 

1 

1 

1 

* 

• 

* 

• 

i 

i 

1 

1 

1 

• 

• 

# 

a 

(V 

c 

© 

c 

r\ 


r' 

r~ 

© 

© 

st 



























© 

♦ 

e 


r\j 

X 

© 


X 

X 

f- 

© 

vC 


tl 
























lT 

© 

© 

4 

■4 

*4 

c 

m 

O' 

m 

c 

< 

4 


o: 
























r- 

© 

O' 


4- 

O 

o 


f- 

jn 

o ru 



























© 

© 

r' 

r- 

IT 

© 


o 

r- 

o 

»>« 

© 

© 


X 









wm 

— * 


_ 


M 

mm 















• 

• 



• 




X 








© 

© 

© 

© 

© 


© 

© 











t 


1 


1 


i 



1 





>- 






1 

I 

1 

1 

1 

1 

1 

1 
























_i 







UJ 

UJ LU 

UJ 

uJ jJ 

L< 

jJ 
























4 


X 

4 

X 

ru 

tr 

c 

4 

© 

rv 

rv 

© 

4 

© 

S 

(V 

a 

4 

© 

ru 


C 















C 

»- 

VJ 

4 

o 

ru 

rv 

iT 

© 

rv 

© 

p^ 

►- 

a 

ru 

© 

s 

rv 

rv 

o 


© 

© 

c* 

n 

© 


lt m 



-4 

o 

if. 


r^- 

O 


© 

X 

o 

O 

© 

r-* 


4 

j-. 

© 


ru 

ru 


© 

© i 

4 


r- 

X 

© 

© 

X 

© 

m 


© 

© 

O 

e 

© 

rw o 

o 

ir in 



© 

4 

_J 

4 

4 

m 

© 

ru 

m 

tP rn 

ru 

rv 




ru 

X 

«n 

4 

4 

4 

4 


f\ 

4 

© 


r* 


e 


r“ 

c 


O 

IT 






*— 


•— 


© 

© 

© 

© 

«£ 

© 


sT 







— 


r-» 

O' 

© 

r— 

mm 

O 

c 

N O IP 

O' in f- 


© 


> 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

r- 

m 

c 

© 

iT 

t-- 

C 

© 

O 



O 


ru 



• 

I 

* 

1 

l 









l 

1 

i 

1 

1 

1 

\ 

l 


i 



1 


i 


1 


l 



£ 






































O 

< 





































a 

X 


X 




































Cl 

UJ 

l- 


IT 

n 

o 

V 

© 

© 

rn 

© 

X 

f*" 

© 

c? 

ru 

o 

'n 

© 


o 

o 

© 

© 

© 

in 

r- 

© 



n. 

ru 

© 

c 


Z 

-J 

UJ 

c 

4 


© 

4 

© 

4 

nn 

iT 

IT 


4 

c 

4 

X 


© 

4 



— 

© 

4 

© 

© 

rv o 

rn 

4 


© 

© 


in 

o 

< 

a 

z 

rv 

m 

ru in 

© 

4 

C 


© 

© 


c 

4 

© 

© 

ru 

© 

ru 

4 


4 

© O 


© 

© IT. 

o 

OJ 

© 


in in o 



D k) O 

o 

r-~ 

r- 

© 

O 


4 

IT 

IT- 

4 


O 

© 

f*. 


o 

© 

C 

o 


O' 



© 

o 


© 

<0 

4 

*> 

© 

r- 

© 

>- 


O 

a. 

© 

r^. 


© 

© 

r*i 

m in 

IT 

if. 

m ir 


© 

© 

p» 

r- 

X 

X 



N © 

© 

e 



© 

© O 

© o 




in 

H- 





















• 

• 














© 


in 






1 

i 

1 

1 

1 

i 

1 

1 










1 



1 


t 


1 


a 



VJ 

© 





































o 

in 







































-j 




































UJ 

cu 

UJ 



no m 

4 

in © 


© 

a 

© 

mm 

ru m 

4 

© 

6 

7 

X 

o 

© 



m ,© 

o 

© «v 

o 

fT 

ru 

r- 

© © 


> 


z 

• 











mm 

— « 


<■— 

mm 

*— 



mm 

rv 

rv 

4 

m 


o in ^ 

4 

© 

£ 


4 

© © 


n 

4 

o 






















fU 


4 

© 

© o 

© 

in © o n 4 

c 


X 

a 

z 






















© 

4 

f- 

o 

© r* 

m 

m 4 

« ru 


ru 


u 
























© IT. 

f- 

m 

© c 

* 

c 

« o 

© © 

c 


4 
























• 

• 

• 

■ 

• 

• 



• 

• 

• 

• 



X 

























m 


<D 

n3 

CL 


153 


Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15.- Continued. 



fi/iM/V Y Z ASHF.U 8t T A YC ZC RG RG/R 




O 

<D 

tr> 

fd 

Cu 


o o 

o => 

co i\i 'vj 

CO CO 

X' X V ru 

. . O' o 

'v "V -c £ 



r\j f\J 

e e 

u r\j ru -w 

0 o x <o 

1 >C »- O o- 

• • »C X 



160 


i 


Figure 15.- Continued. 
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Figure 15.- Continued. 
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Figure 15.- Concluded. 










Figure 17.- Measured and predicted circumferential 
pressure distribution on an ogive-cylinder body 
at Mco = 1.6, a c = 20°. 
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Figure 21.- Measured and predicted downwash distribution 
on the lee side of an ogive cylinder 
at M„ = 2, a _ = 15°, x/D = 10 
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Figure 25.- Measured and predicted circumferential 
distribution on a 3:1 elliptic cross section 
missle at Moo - 2.5, a c - 20°, x/L = 0.60. 
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Figure 26.- Concluded. 
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EXPERIMENT (Ref. 18) 
Theory 



Figure 32.- Measured and predicted pressure distribution on a lobed cross 

section body at = 1.7, a c = 20.36°. 
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problem for solution. Comparison of measured and predicted 
surface pressure distributions, flow field surveys, and aero- 
dynamic chararacteristics are presented for bodies with 
circular and noncircular cross-sectional shapes. 
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